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SECTION  I 
INTRODUCTION 


A.  OBJECTIVE 

The  objective  of  this  program  was  to  determine  the  effects  of  VOC  con¬ 
centration,  temperature,  species  competition,  and  humidity  on  the  sorption 
performance  (capacity  and  kinetics)  of  activated  carbons  and  possibly  other 
sorbents  for  the  removal  of  selected  VOCs.  VOCs  of  interest  include  ben¬ 
zene,  ethylbenzene,  toluene,  xylenes,  trichloroethylene,  and  1,1-dichloro- 
ethylene.  Concentrations  of  interest  are  in  the  range  of  30  to  300  mg/m3 
(nominally  7  to  70  ppmv). 

B.  BACKGROUND 

Air  stripping  is  a  widely  used  cleanup  technique  for  removing  volatile 
organic  contaminants  (VOCs)  from  groundwater.  The  effluent  air  stream  is 
usually  vented  directly  to  the  atmosphere.  Although  the  resulting  VOC  con¬ 
centrations  in  the  air  are  relatively  low,  there  is  increasing  pressure 
from  both  state  and  federal  agencies  to  reduce  the  air  pollution  potential 
of  such  processes.  The  Air  Force  Is  considering  the  use  of  such  stripping 
processes  to  remove  VOCs  such  as  benzene,  ethylbenzene,  and  trichloro¬ 
ethylene  from  groundwater. 

Adsorption  of  a  contaminant  from  a  stream  of  gas  passing  through 
either  a  fixed  or  fluidized  bed  of  sorbent  is  a  dynamic  process  that 
depends  on: 

1.  The  equilibrium  distribution  of  sorbate  in  the  vapor  and  adsorbed 
phases. 

2.  The  dynamics  of  vapor  transport: 

(a)  to  particle  surfaces,  and 

(b)  within  the  micropore  structure  of  thw  sorbent  particles. 

3.  The  rate  of  flow  through  the  bed. 
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The  equilibrium  between  sorbate  and  the  adsorbed  phase  Is  a  function 
of  the  concentration  of  the  sorbate  In  the  gas  phase;  this  Is  usually 
expressed  In  terms  of  an  adsorption  Isotherm.  Adsorption  Isotherms  can 
have  several  different  shape  characteristics  depending  on  the  specific 
types  of  Interactions  (surface  adsorption,  multilayer  adsorption,  pore  con¬ 
densation,  etc.),  and  these  can  generally  be  categorized  In  terms  of 
Brunauer's  classification  system  (Reference  1).  For  example,  the  adsorp¬ 
tion  of  water  on  activated  carbons  usually  follows  Brunauer's  (BET)  Type  II 
isotherm  with  hysteresis,  illustrated  In  Figure  1  In  the  form  of  a  plot  for 
the  amount  adsorbed  (arbitrary  units)  versus  the  ratio  of  the  partial  pres¬ 
sure  to  the  equilibrium  vapor  pressure.  This  type  of  Isotherm  shows  three 
distinct  regions: 


1.  At  low  relative  pressures  (below  the  knee),  adsorption  occurs 
primarily  on  the  free  surfaces  of  the  sorbent. 

2.  At  Intermediate  pressures,  the  Isotherm  Is  relatively  flat,  with 
only  multilayer  adsorption  occurring. 

3.  At  some  critical  relative  pressure  (0.6  In  Figure  1),  depending  on 
the  micropore  size  distribution,  capillary  condensation  begins  to 
occur,  with  the  result  that  the  micropores  become. filled  with 
liquid  water.  This  process  Is  possible  because  of  the  small  size 
of  the  pores  and  hence  the  small  radius  of  curvature  of  the  liquid 
surface. 


This  behavior  may  inhibit  the  kinetics  of  adsorption  of  other  com¬ 
ponents  (especially  those  present  at  low  concentrations)  because  of  the 
fact  that  pore  surfaces  can  become  blocked  by  liquid  water.  In  such  cases, 
adsorption  of  the  trace  components  may  still  occur  if  they  are  soluble  in 
water,  but  this  adds  more  resistance  to  the  overall  transport  process. 

Transport  processes  within  the  sorbent  bed  are  determined  primarily  by 
the  particle  size  and  the  pore  structure  of  the  particles.  Large  particles 
yield  large  Intraparticle  spaces,  and  decrease  the  efficiency  of  contact 
between  the  gas  and  the  sorbent  bed.  On  the  other  hand,  very  small  Intra- 
particle  spaces  cause  high  pressure  drops  across  the  bed,  and  some  trade¬ 
offs  are  necessary.  Transport  within  a  particle  occurs  solely  by  diffu¬ 
sion.  Thus,  large  pores  provide  the  best  access  to  the  internal  surfaces 
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Figure  1.  Typical  BET  Type  II  Isotherm. 


of  the  particles.  On  the  other  hand,  smaller  pores  provide  more  surface 
area  per  unit  porosity  and  therefore  provide  greater  sorption  capacity 
per  unit  weight  of  sorbent.  With  pore  sizes  that  are  smaller  than  about 
twice  the  mean  free  path  of  the  sorbate  molecules  however,  the  effective 
diffusivitles  of  the  sorbate  molecules  become  smaller  than  their  free  space 
diffusivities,  and  their  access  to  the  sorbent  surface  becomes  kinetically 
limited.  In  the  extreme  examples  of  such  cases,  the  sorbent  may  have  a 
large  equilibrium  capacity  for  the  sorbate,  but  will  not  perform  well  under 
dynamic  sorption  conditions,  i.e.,  when  the  flow  rate  through  the  bed  is 
faster  than  the  speed  at  which  the  sorbent  can  accept  the  sorbate. 

In  practical  testing  and  evaluation  of  sorbent  beds,  it  is  therefore 
important  that  both  the  equilibrium  and  kinetic  parameters  be  evaluated  so 
that  appropriate  tradeoffs  can  be  considered  for  development  of  a  system 
that  is  appropriate  for  a  given  task.  The  most  commonly  used  engineering 
approach  involves  the  use  of  the  Wheeler  model  (or  one  of  its  several  modi¬ 
fications)  of  bed  permeation  (Reference  2).  Experimentally,  the  approach 
involves  the  measurement  of  bed  breakthrough  times  (tb)  for  several  differ¬ 
ent  bed  weights  (W),  and  the  data  are  then  correlated  using 


where  We  is  the  equilibrium  bed  loading  (g/g)  at  the  inlet  concentration  Ci 
(g/L),  Q  is  the  volumetric  flow  rate  (L/min),  d  is  the  bed  density  (g/L),  k 
is  the  overall  rate  constant  (min-1),  and  C  is  the  outlet  gas  concentration 
at  time  tb-  We/C-j  defines  a  point  on  the  adsorption  isotherm  at  C-j,  and  at 
very  low  concentrations,  where  the  isotherm  is  approximately  linear,  this 
ratio  defines  the  slope  of  the  isotherm.  Inasmuch  as  this  slope  is  not 
constant  over  the  whole  isotherm  (see  Figure  1),  it  is  important  that  bed 
performance  be  measured  at  concentrations  expected  to  prevail  under  normal 
operating  conditions.  In  the  case  of  VOC  removal  from  humid  air  streams, 
it  is  therefore  necessary  to  determine  sorption  behavior  at  concentrations 
in  the  range  of  interest  (30  to  300  mg/m^).  For  the  VOCs  of  interest  to 
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this  program,  this  concentration  range  will  yield  relative  pressures,  P/P0 
of  the  order  of  10"8  to  10-5. 

The  BET  isotherm,  illustrated  above  for  the  adsorption  of  water  vapor 
on  carbon,  is  not  well  suited  for  description  of  the  adsorption  process  for 
VOCs  at  these  low  levels  of  P/P0.  For  such  cases,  the  Dubinin-Radushkevich 
isotherm,  viz. 

i .  .  [.  Sf  (,*•  y] 

has  proven  of  much  greater  utility  (Reference  3).  In  this  equation,  wv  is 
the  volume  of  condensed  adsorbate,  W0  is  maximum  volume  available  for  sor- 
bate  condensation,  B  is  a  constant  that  is  characteristic  of  the  adsorbent, 
and  b  is  an  affinity  coefficient  that  compares  the  strength  of  the  adsorp¬ 
tive  Interaction  with  some  reference  adsorbate.  Jonas  et  al  (op.  cit.) 
have  shown  that  values  of  b  for  different  sorbates  can  be  correlated  with 
their  electronic  polarizabilities.  Therefore,  once  a  particular  adsorbent 
has  been  characterized  with  respect  to  a  reference  adsorbate  (i.e.,  once  B 
is  determined  for  the  adsorbent),  then  the  adsorption  behavior  of  that  sor¬ 
bent  towards  other  sorbates  can  be  predicted.  More  recently,  Jonas  and 
co-workers  (Reference  4)  have  extended  this  approach  to  prediction  of  car¬ 
bon  adsorption  performance  for  binary  vapor  mixtures. 

C.  SCOPE 

Ten  commercially  available  carbon  sorbents  were  evaluated  for  ability 
to  adsorb  selected  VOCs.  The  more  promising  of  these  sorbents  were  sub¬ 
jected  to  detailed  adsorption  testing  as  a  function  of  temperature,  flow 
rate,  sorbate  composition,  and  relative  humidity.  Results  were  correlated 
using  the  Oubinin-Radushkevich  isotherm  and  the  Wheeler  bed  permeation 
model,  and  a  computer-based  predictive  model  was  developed.  Based  on 
results  obtained,  an  economic  assessment  was  made  for  several  possible 
engineering  options  for  carbon-based  cleanup  of  VOC  contaminated  ground- 
water.  It  is  concluded  that  air  stripping  followed  by  purification  of  the 
stripper  air  Is  more  economical  than  direct  contact  of  the  well  water  with 
activated  carbon. 
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SECTION  II 

EXPERIMENTAL  APPROACH 


A.  GENERAL 


The  approach  taken  in  this  program  can  be  divided  Into  a  series  of 
subtasks: 


•  -  Adsorbent  screening  and  selection.  Ten  activated  carbons  were 

considered  initially  as  candidate  sorbents.  These  sorbents  were 
subjected  to  a  series  of  measurements  to  determine  their  relative 
adsorption  capacities,  adsorption  kinetics,  and  water  adsorption 
Isotherms.  Based  on  the  results  of  these  tests  and  considerations 
of  sorbent  costs,  the  list  of  sorbents  to  be  considered  further 
was  limited  to  the  five  that  appeared  to  be  most  promising.  This 
list  was  further  reduced  during  subsequent  tests  until  only  one 
sorbent  remained  for  the  final  series  of  experiments. 

•  -  Single-component  studies.  Adsorption  studies  were  conducted  using 

representative  VOCs  one  at  a  time.  These  studies  were  conducted 
using  various  concentrations  at  ambient  temperature  and  low 
(<5  percent)  humidity  to  provide  an  optimum  performance  baseline 
set  of  data. 

•  -  Multiple-component  studies.  The  adsorption  behavior  of  various 

combinations  of  the  VOCs  of  interest  were  determined  at  low 
humidities. 

•  -  Humidity  effect  studies.  The  adsorption  of  VOCs  of  interest  was 

studied  at  various  humidity  levels.  These  studies  were  conducted 
with  both  single-  and  multiple-component  feedstocks  at  ambient 
temperature.  Finally  the  effect  of  increased  temperature  on 
adsorption  of  multiple-component  mixtures  of  VOCs  from  humidified 
feedstocks  was  evaluated. 

•  -  Economic  assessment.  A  preliminary  process  design  and  economic 

analysis  was  performed  to  obtain  capital  and  operating  costs  for 
treatment  of  well  water  contaminated  with  VOCs  by  activated  carbon 
adsorption.  A  total  of  nine  treatment  options  were  Included  in 
the  analysis  to  assess  vapor-phase  adsorption  versus  liquid-phase 
adsorption;  effects  of  contaminant  levels  in  the  well  water;  and 
effects  of  relative  humidity  levels  in  the  air  stream  feed  to  the 
absorber. 


Throughout  the  course  of  the  program,  data  were  analyzed  In  terms  of 
the  Dubinin-Radushkevich  and  Wheeler  (DRW)  models.  To  facilitate  these 


analyses,  a  computer  model  incorporating  the  DRW  concepts  and  using  ideal 
solution  theory  to  estimate  multiple-component  effects  was  developed.  The 
mathematical  development  of  this  model  is  described  in  Appendix  A. 

B.  APPARATUS 

A  schematic  representation  of  the  experimental  apparatus  used  for  most 
of  the  experiments  conducted  on  this  program  is  shown  in  Figure  2.  With 
this  system,  a  metered  flow  of  gas  containing  the  VOCs  is  passed  through 
the  sorbent  bed  at  controlled  levels  of  humidity  and  temperature.  The 
effluent  then  passes  through  a  gas  sampling  valve  to  a  continuously  oper¬ 
ating  flame  ionization  detector  (FID).  With  this  configuration,  break¬ 
through  of  the  VOCs  was  detected  by  the  FID  detector.  When  using  multiple- 
component  sorbate  mixtures,  repeated  samples  of  the  effluent  stream  were 
automatically  routed  to  a  gas  chromatograph  for  analysis.  Calibrations  of 
the  FID  and  the  gas  chromatograph  were  performed  by  bypassing  the  sorbent 
bed.  With  this  system,  complete  breakthrough  curves  could  be  generated, 
and,  when  necessary,  desorption  curves  could  be  observed.  For  those 
experiments  requiring  high  humidity,  a  Nafion  dryer  was  installed  at  the 
bed  outlet.  This  avoided  problems  associated  with  condensation  and  subse¬ 
quent  aqueous  phase  absorption  of  VOCs  and  interference  of  water  with  the 
analysis  procedures. 

C.  SORBENT  BEDS 

Sorbent  beds  were  prepared  in  I/4-inch  stainless  steel  tubes  con¬ 
taining  0. 1-0.8  grams  of  carbon  depending  on  the  density  of  the  carbon  and 
the  depth  of  bed  desired.  Before  each  use,  beds  were  baked  at  180- 
200  degrees  Centigrade  with  nitrogen  passing  through  the  bed.  Baking  was 
continued  until  no  evidence  of  hydrocarbon  elution  was  seen  with  the  FID. 

In  initial  tests,  fresh  beds  were  prepared  for  each  experiment.  These 
tests  Included  thermal  desorption  measurements  made  at  the  completion  of 
each  adsorption  run.  Subsequent  tests  with  the  thermally  regenerated  beds 
showed  no  detectable  differences  in  adsorption  performance.  A  limited 
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number  of  regeneration  trials  were  conducted  using  steam  regeneration  of 
the  beds,  but  no  particular  advantage,  other  than  the  use  of  lower  temper¬ 
atures  was  found  for  this  procedure.  Subsequent  experiments  were  conducted 
using  thermally  regenerated  beds,  with  occasional  comparisons  being  made 
with  fresh  beds.  Deterioration  of  bed  performance  occurred  only  once, 
after  several  tens  of  regenerations. 

D.  CARBONS 

Activated  carbons  considered  Initially  for  this  program  are  cited  in 
Table  1.  These  Included  both  conventional  carbons  and  recently  developed 
specialty  carbons. 

TABLE  1.  LIST  OF  SORBENTS  USED  IN  SCREENING  TESTS 

Sorbent  Type  Comments 

Activated  Carbons: 

SK-4 
LCL 
WV-B 
CT 

Specialty  Carbons: 

Graphpac 
Carbosphere 

Spherocarb 
Carbopack 
Carbosleve  G 
Carbosleve  S-ll 

E.  ADSORPTION  ISOTHERMS 

Low-concentration  data  for  the  Oubinin-Radushkevich  isotherms  were 
derived  directly  from  the  bed-permeation  experiments  described  above. 


graphltlzed  carbon 
carbon  "mole  sieve" 

carbon  "mole  sieve" 
graphltlzed  carbon 
low-density  carbon 
high-density  carbon 


area  -  100  m2/g 
1000  n»2/g;  low  water 
retention 

similar  to  Carbosphere 
similar  to  Graphpac 

ca.lOOOm^/g 


nut  shell 
coconut  shell 
coal  base 
coconut  shell 


* 
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However,  high-concentration  (near  saturation)  capacities  for  the  VOCs  and 
water-adsorption  Isotherms  were  more  conveniently  determined  gravlmetrl- 
cally.  For  these  measurements,  a  Cahn  Model  RG  recording  microbalance  was 
used  with  the  following  procedure: 


1.  A  sample  of  the  carbon  was  placed  on  the  balance  pan,  and  was 
outgassed  at  about  150  degrees  Centigrade  in  a  stream  of  dry 
nitrogen. 

2.  After  cooling  the  sample  to  25  degrees  Centigrade,  the  gas  flow 
to  the  sample  area  was  switched  to  either  a  humidified  stream  of 
nitrogen  or  nitrogen  containing  a  high  concentration  of  the  VOC  of 
interest. 

3.  The  sample  was  allowed  to  equilibrate  with  the  gas  stream,  as 
judged  by  the  balance  output.  During  this  time,  the  sorbate  con¬ 
centration  was  measured  using  either  a  relative  humidity  meter  for 
water  concentrations  or  a  gas  chromatograph  to  monitor  the  VOCs. 
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4.  Steps  2-3  were  then  repeated  with  successively  higher  sorbate 
concentrations. 


SECTION  III 

RESULTS  AND  DISCUSSION 


A.  SCREENING  TESTS 

Each  of  the  10  carbons  cited  in  Table  1  was  subjected  to  an  initial 
set  of  screening  tests  that  consisted  of  determining:  (a)  breakthrough 
curves  at  low  (<5  percent)  and  high  (85-90  percent)  humidities  for  both 
1,1-dlchloroethylene  and  benzene;  (b)  thermal  desorption  efficiencies;  and, 
(c)  water  adsorption  isotherms  at  ambient  temperature.  Detailed  results  of 
these  measurements  are  given  in  Appendix  B.  Table  2  summarizes  these  mea¬ 
surements  in  terms  of:  (a)  relative  adsorption  capacities;  (b)  relative 
standard  deviations  for  the  breakthrough  curves;  (c)  water  adsorption 
capacity  at  saturation;  and  (d)  relative  humidity  required  for  50  percent 
saturation  of  the  carbon.  Capacities  were  calculated  both  by  total  inte¬ 
gration  of  the  desorption  curves  and  by  using  the  50  percent  breakthrough 
time  and  assuming  a  symmetrical  breakthrough  curve.  The  former  approach  is 
considered  more  accurate,  but  the  two  methods  agreed  to  within  a  few  per¬ 
cent  for  all  runs. 

The  data  in  Table  2  Indicate  that  the  two  Carbosieve  carbons  consis¬ 
tently  yield  the  highest  capacities,  regardless  of  the  exposure  condition. 
The  kinetics  for  these  two  carbons,  as  judged  from  the  standard  deviations, 
also  appear  reasonably  good,  although  there  is  more  scatter  apparent  in 
these  numbers  than  in  the  capacity  data.  CT  was  clearly  the  best  of  the 
nonspecialty  carbons,  while  WV-B  was  one  of  the  poorer  performers.  The 
capacities  of  Graphpac  GB  and  Carbopack  C  were  so  low  that  use  of  these  two 
carbons  was  rejected  very  early.  Therefore,  not  all  of  the  screening  tests 
were  completed  with  these  carbons. 

The  water-adsorption  Isotherms  for  the  eight  more  active  carbons  (see 
Appendix  B)  are  generally  similar,  but  some  differences  can  be  noted.  For 
example,  WV-B  absorbs  more  water  at  70  percent  humidity  than  does  Sphero- 
carb  at  near  saturation.  CT,  SK-4,  Spherocarb,  Carbosphere,  and  Carbosieve 
S-ll  all  have  water  capacities  that  are  significantly  lower  than  those  of 
Carbosieve  G,  LCL,  and  WV-B. 
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TABLE  2.  RELATIVE  SORBENT  PERFORMANCE 


Sorbent 

C/Cmaxa 

smin/sl> 

RH 

Cwater*  50% 

mg/g  Capacity 

OCE 

Benzene 

DCE 

Benzene 

Wo 

D 

W 

D 

W 

D 

W 

SK-4 

0.63 

0.58 

0.57 

0.63 

0.41 

0.32 

0.32 

0.18 

306 

49 

Carbosieve 

S-ll 

1.00 

1.00 

0.98 

1.00 

1.00 

0.65 

0.20 

0.38 

346 

38 

Graphpac  GB 

0 

0 

0 

0 

NA 

NA 

NA 

NA 

NA 

NA 

Carbosphere 

0.68 

0.70 

0.77 

0.69 

0.35 

0.34 

0.50 

0.64 

325 

55 

Carbopack  C 

0 

0 

0 

0 

NA 

NA 

NA 

NA 

NA 

NA 

WV-B 

0.26 

0.38 

0.58 

0.48 

0.38 

0.31 

0.16 

1.00 

528 

70 

Spherocarb 

0.63 

0.85 

0.76 

0.70 

0.48 

0.30 

1.00 

0.12 

255 

60 

LCL 

0.53 

0.45 

0.71 

0.62 

0.15 

1.00 

0.32 

0.72 

410 

62 

CT 

0.76 

0.94 

0.73 

0.73 

0.25 

0.41 

0.41 

0.51 

323 

56 

Carbosieve  G 

0.76 

0.75 

1.00 

0.87 

0.58 

0.50 

0.52 

0.82 

476 

44 

aC  =  sorption  capacity  under  test  conditions;  Cmax  =  maximum  capacity 
observed  in  this  series. 


bs  =  standard  deviation  of  breakthrough  curve;  smin  =  minimum  standard 
deviation  observed. 

CD  =  low  humidity  runs. 

dW  =  high  humidity  runs. 

Figure  3  shows  a  scattergram  combining  all  of  the  data  for  the  desorp¬ 
tion  efficiencies  measured  during  the  screening  tests.  No  consistent  dif¬ 
ferences  in  desorption  efficiency  were  noted  between  individual  carbons. 
Points  circled  in  this  figure  were  deleted  from  the  numerical  analysis  as 
being  more  than  two  standard  deviations  from  the  mean.  The  two  low 
rejected  points  were  observed  with  carbons  LCL  and  SK-4,  and  were  not 
repeated  in  other  runs  with  these  sorbents.  The  three  high  rejected  points 
occurred  with  benzene  desorption  from  Spherocarb  and  Carbosieve  S-ll.  In 
the  case  of  the  former  sorbent,  benzene  is  formed  from  decomposition  of 
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Figure  3.  Thermal  Desorption  Efficiencies. 


incompletely  pyrolyzed  polymer.  The  high  desorption  efficiency  with 
Spherocarb  Is,  therefore,  not  an  unexpected  result. 

Five  of  the  candidate  carbons  were  subjected  to  both  thermal  desorp¬ 
tion  and  multiple  steam  regeneration  cycles  to  determine  If  such  processes 
would  adversely  affect  their  performance.  Results  of  these  measurements 
are  shown  In  Table  3.  This  series  of  treatments  had  little  effect  on  the 
performance  of  any  of  the  carbons.  Indeed,  the  observed  changes  are,  for 
the  most  part,  well  within  the  expected  experimental  error  of  the  measure¬ 
ments.  Consequently,  most  of  the  remaining  experiments  were  conducted  with 
repeated  usage  of  the  same  sorbent  cartridges  rather  than  with  freshly  pre¬ 
pared  cartridges  for  each  experiment.  Spot  checks  on  the  reliability  of 
this  procedure  were  conducted  throughout  the  remainder  of  the  program,  and 
the  frequency  of  such  repeat  runs  was  greatest  during  the  early  stages  of 
successive  tasks. 

TABLE  3.  EFFECT  OF  THERMAL  DESORPTION  AND  STEAM  TREATMENT 
ON  SORBENT  PERFORMANCE  (BENZENE  AT  300  mg/m^) 


Sorbent 

Capacity,  g/ga 
(virgin) 

Capacity,  g/g 
(treated)15 

Percentage 

Change 

SK-4 

0.114 

0.117 

+2.6 

Carbosieve  S-ll 

0.195 

0.180 

-7.7 

Spherocarb 

0.149 

0.151 

+1.3 

CT 

0.143 

0.141 

-1.4 

Carbosieve  G 

0.200 

0.191 

-4.5 

Mean 

-1.9  +/-4.2 

Capacity  based  on  50  percent  breakthrough  time  and  Input 
concentration  level.  (300  mg/nn  x  flow  rate  x  tso/bed  weight.) 

bV1rg1n  samples  exposed;  thermally  desorbed  at  160  degrees  Centigrade; 
steam  desorbed  at  100-140  degrees  Centigrade;  and  reexposed  to  benzene. 

Based  on  the  results  of  these  screening  tests,  five  carbons,  CT, 
Carbosieve  S- 11,  SK-4,  Carbosleve  G,  and  Spherocarb,  were  selected  for 
further  evaluation. 


B.  REFERENCE  DATA  FOR  MODEL  DEVELOPMENT 


The  five  carbons  cited  above  were  used  to  develop  a  reference  data  set 
for  use  with  the  DRW  model.  In  these  experiments,  adsorption  breakthrough 
times  (defined  as  the  time  when  the  eluent  concentration  rises  to  10  per¬ 
cent  of  the  input  concentration)  and  adsorption  capacities  were  measured  as 
a  function  of  bed  depth,  feed  rate,  and  input  concentration  using 

1.1- dichloroethylene  as  the  reference  sorbate.  As  noted  in  the  experimen¬ 
tal  section,  additional  capacity  measurements  were  made  at  very  high 
(P/Po  -  1)  1,1-dichloroethylene  concentrations  using  the  gravimetric  proce¬ 
dure.  The  combined  data  were  then  used  to  derive  appropriate  DR  parameters 
for  each  of  the  five  carbons.  A  detailed  tabulation  of  the  results 
obtained  in  this  series  of  experiments  is  given  in  Appendix  C.  Initial 
estimates  of  the  rate  constants  indicated  considerable  uncertainty  in  that 
parameter  (see  Appendix  C).  Additional  measurements  of  bed  depth  effects 
were,  therefore.  Incorporated  Into  the  concentration  effect  study  In  an 
effort  to  gain  a  better  overall  representation  of  the  reference  case  with 

1.1- dichloroethylene.  The  apparent  rate  constants  were  quite  large  com¬ 
pared  with  those  cited  by  Jonas  et  al.  This  can  be  attributed  to  the  much 
smaller  carbon  particles  being  used  in  the  current  study.  Jonas  found 
apparent  rate  constants  of  about  1000  min-*  with  carbon  particles  having 
mean  diameters  of  0.268  cm.  The  current  work  is  being  conducted  with 
60-80  mesh  (d  =  0.0216  cm)  particles.  Because  the  transport  process  to  and 
within  the  particles  Is  diffusion-limited,  the  rate  of  transport  should  be 
inversely  proportional  to  the  particle  diameter.  Therefore,  we  would 
expect  the  rate  constants  for  the  current  work  to  be  of  the  order  of 
10,000  m1n~l,  with  minor  variations  due  to  the  internal  structures  of  the 
different  carbons.  To  make  full  use  of  all  the  data,  the  Wheeler  equation 
was  used  along  with  the  measured  capacities  for  each  experiment  to  calcu¬ 
late  an  apparent  rate  constant  for  each  experiment.  These  values  were  then 
grouped  for  each  sorbent,  and  the  best  value  was  determined  for  each 
sorbent.  These  results  are  shown  in  Table  4. 
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TABLE  4.  SUMMARY  OF  RATE  CONSTANTS 


! 

i 


Sorbent 

Rate  Constant,  min"l 

RSD,  % 

Carbosieve  S-ll 

2.14  x  104 

21 

CT 

1.86  x  104 

22 

SK-4 

1.23  x  104 

19 

Carbosieve  G 

1.70  x  104 

18 

Spherocarb 

1.20  x  104 

26 

The  apparent  rate  constants  are  of  the  right  magnitude  compared  with 
those  measured  by  Jonas.  Although  the  uncertainties  are  large,  with  this 
large  rate  constant,  the  contribution  of  the  kinetic  term  to  the  break¬ 
through  time  is  relatively  small,  i.e.,  the  breakthrough  time  is  rela¬ 
tively  insensitive  to  the  exact  value  of  the  rate  constant.  (Conversely, 
the  calculated  rate  constants  are  sensitive  to  minor  variations  in  the 
breakthrough  times.) 

Figure  4  shows  an  example  of  the  plot  used  to  derive  the  DR  para¬ 
meters  for  CT.  Table  5  summarizes  the  DR  parameters  derived  from  the  data 
included  in  this  part  of  the  study.  Table  6  summarizes  the  precision  and 
accuracy  of  application  of  the  DR  isotherm  to  the  reference  data  for 
1,1-dichloroethylene. 

TABLE  5.  DUBININ-RADUSHKEVICH  PARAMETERS 


Sorbent 

In  Wo 

+/-a 

k  x  108 

+/-  (x  108)a 

Carbosieve  S-II 

6.26 

0.039 

5.03 

0.10 

CT 

6.05 

0.056 

5.02 

0.16 

SK-4 

5.98 

0.055 

5.60 

0.16 

Carbosieve  G 

6.36 

0.045 

5.92 

0.12 

Spherocarb 

6.08 

0.049 

5.39 

0.14 

aProbable  error 


Figure  4.  Dubinin-Radushkevich  Plot  for  1,1-Dichloroethylene  on  CT 


TABLE  6.  SUMMARY  OF  PRECISION  AND  BIAS  OF  DRW  REFERENCE  DATA 


Sorbent 

*10 

Capacity, 

mq/q 

Pa 

Bb 

P 

B 

Carbosieve  S-ll 

9.31 

-1.49 

10.0 

-2.56 

CT 

19.3 

5.41 

11.5 

-2.27 

SK-4 

27.6 

8.92 

9.37 

-0.16 

Carbosieve  G 

14.9 

10.1 

7.80 

0.54 

Spherocarb 

24.3 

-11.4 

10.2 

-0.95 

“Precision  defined  as  root  mean  square  percentage  dif¬ 
ference  between  predicted  and  measured  tio  or  capacity. 


bB1as  defined  as  average  percentage  difference  between 
predicted  and  measured  t 10  or  capacity.  (A  negative  bias 
Indicates  low  predicted  values.) 


C.  SINGLE  COMPONENT  VOC  ADSORPTION 

Additional  bed  permeation  experiments  were  conducted  using  various 
combinations  of  the  five  candidate  sorbents  at  various  bed  depths  for 
adsorption  of  benzene,  trichloroethylene,  toluene,  o-xylene,  and  ethyl¬ 
benzene.  Detailed  results  for  these  experiments  are  given  in  Appendix  D. 
The  DRW  parameters  derived  from  the  reference  data  set  with  1,1-dlchloro- 
ethylene  were  used  to  estimate  values  of  t^g  and  the  adsorption  capacities 
for  each  run.  The  calculated  and  measured  results  were  then  compared. 
Application  of  the  DRW  model  to  compounds  other  than  the  reference  chemical 
is  described  In  Appendix  A.  Briefly,  two  adjustments  to  the  reference 
parameters  are  required: 

1.  The  saturation  capacity.  Wo,  must  be  adjusted  to  accommodate  the 
difference  in  density  between  the  reference  chemical  and  the  test 
chemical . 

2.  The  affinity  coefficient  for  the  test  chemical  must  be  calculated 
relative  to  that  of  the  reference  chemical  using  the  ratios  of 
their  electronic  polarizabil ities. 
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Appropriate  parameters  for  estimation  of  the  affinity  coefficients  of 
the  test  VOCs  are  given  in  Table  7.  The  results  of  application  of  the  DRW 
model  to  the  single  component  data  are  summarized  in  Table  8. 


TABLE  7.  S0R8ATE  PROPERTIES  FOR  DR  CORRELATION 


Chemical 

Mwa 

P0  mg/m3 
(x  10-6) 

@  25C 

db 

nc 

ped 

pe/peDCE 

1,1-Dichloro- 

96.94 

31.9 

1.218 

1.425 

20.35 

1.00 

ethylene 

Trichloro- 

131.39 

5.07 

1.464 

1.477 

25.37 

1.25 

ethylene 

o-Xylene 

106.17 

0.353 

0.880 

1.506 

35.81 

1.76 

Toluene 

92.15 

1.34 

0.867 

1.496 

31.06 

1.53 

Ethylbenzene 

106.17 

0.546 

0.867 

1.496 

35.77 

1.76 

Benzene 

78.12 

3.82 

0.879 

1.501 

26.20 

1.29 

aMolecular  weight 
bLiqu1d  density,  g/cc 
cRefract1ve  index  of  liquid 
^Electronic  polarizability 


Comparing  the  data  of  Table  6  with  those  in  Table  8  shows  that  the 
indicated  precision  are  probably  inherent  in  the  experimental  technique. 
This  could  be  due  to  variations  in  bed  packing  densities.  In  later  experi¬ 
ments  where  the  number  of  variables  (principally  the  number  of  different 
carbon  beds)  was  much  more  limited,  the  precision  and  apparent  bias  were 
considerably  improved  over  the  data  shown  in  Table  8.  However,  even  the 
range  of  bias  Indicated  by  Table  8  would  probably  not  be  unacceptable  for 
purposes  of  engineering  estimates  of  bed  capacities.  The  biases  listed  in 
Table  8  for  the  tjo  estimates  parallel  those  for  the  capacity  estimates. 
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TABLE  8.  APPLICABILITY  OF  DRW  MODEL  TO  TEST  VOCs 


Sorbent 

tio 

Capacity, 

mg/9 

pa 

Bb 

P 

B 

Benzene 

12.2 

7.87 

9.87 

6.80 

Trichloroethylene 

29.9 

14.0 

22.2 

13.2 

o-Xylene 

31.8 

-26.8 

31.1 

-26.6 

Toluene 

21.8 

-17.4 

18.1 

-14.2 

Ethylbenzene 

27.5 

-22.7 

24.3 

-19.3 

Precision  defined  as  root  mean  square  percentage  dif¬ 
ference  between  predicted  and  measured  tio  or  capacity. 


bBias  defined  as  average  percentage  difference  between 
predicted  and  measured  tio  or  capacity.  (A  negative  bias 
indicates  low  predicted  values.) 

This  suggests  that  the  source  of  the  bias  is  predominantly  in  the  corre¬ 
lation  of  the  affinity  coefficients  with  electronic  polarizabilities,  and 
that  a  reexamination  of  the  parameters  used  in  this  correlation  might  lead 
to  improved  overall  prediction  of  capacities. 

D.  MULTIPLE-COMPONENT  STUDIES 

While  considerable  work  has  been  done  to  characterize  adsorption  of 
individual  sorbates  on  a  wide  variety  of  sorbents,  relatively  little  effort 
has  been  directed  towards  the  simultaneous  adsorption  of  multiple  sorbates. 
Historically,  the  adsorption  capacities  of  carbons  and  other  sorbents  have 
been  characterized  in  terms  of  one  or  more  individual  sorbates.  Such  data 
are  then  used  as  "best-case"  examples  for  estimating  bed  performance. 


When  more  than  one  sorbate  is  present  in  the  feedstock,  occurrences  of 
both  competitive  and  cooperative  types  of  phenomena  can  make  the  overall 
description  of  the  adsorption  process  quite  complicated.  Consider  two 
extreme  cases  Involving  only  two  sorbates: 

Case  1.  Two  completely  immiscible  sorbates  --  Because  these  sorbates 
are  completely  immiscible  they  cannot  condense  on  the  same 
volume  of  sorbent  at  the  same  time.  They  must, therefore , 
occupy  separate  zones  of  the  sorbent  bed,  with  the  leading 
2 one  being  occupied  by  the  sorbate  less  strongly  held. 
Furthermore,  continued  addition  of  the  feedstock  to  the  bed 
will  result  in  expansion  of  the  trailing  zone  and  concomitant 
displacement  of  the  leading  zone.  As  the  feed  continues  to 
the  leading  zone,  its  movement  will  be  accelerated.  This 
type  of  behavior  Is  approximated  by  combinations  of  water  and 
many  hydrocarbons,  and  the  phenomenon  is  the  basis  for  steam 
desorption  of  carbons. 

Case  2.  Two  completely  miscible  sorbates  —  If  two  sorbates  are  com¬ 
pletely  miscible  and  their  adsorption  potentials  are  identi¬ 
cal  (i.e.,  they  have  the  same  chemical  functionality  and 
vapor  pressures),  they  will  coadsorb  on  the  same  volume  of 
sorbent  without  separating  into  zones.  An  example  of  this 
type  of  behavior  is  seen  with  the  meta-  and  parai somers  of 
xylene,  which  are  extremely  difficult  to  separate  by  chroma¬ 
tographic  means. 

With  real  sorbate  systems,  even  when  water  is  present  in  the  feed¬ 
stock,  there  Is  usually  some  mutual  solubility  of  the  sorbates,  resulting 
in  zone  formation  and  at  least  some  limited  coadsorption.  Even  with 
chemically  similar  sorbates,  such  as  benzene  and  toluene,  separation  into 
zones  and  accelerated  movement  of  the  leading  zone  are  expected  because  of 
differences  In  vapor  pressure.  This  implies  that  the  greater  the  number  of 
sorbates,  the  greater  will  be  the  differences  between  retention  of  indi¬ 
vidual  sorbates  in  the  mixture  and  their  retention  as  single  components. 


I 

i 
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In  other  words.  In  more  complex  systems,  less  meaning  Is  given  to  estimates 
based  on  single  component  reference  systems.  The  only  useful  approach  Is 
to  employ  the  OR  isotherm  to  estimate  adsorption  capacities  In  each  zone, 
using  Ideal  solution  theory  to  calculate  the  adsorbate  compositions.  The 
details  of  this  approach  are  described  In  Appendix  A. 

Experiments  in  this  portion  of  the  program  involved  a  single  carbon, 
CT,  and  various  mixtures  of  the  six  test  VOCs.  The  target  matrix  of  tests, 
based  on  representative  conditions  at  selected  Air  Force  Bases,  is  given  In 
Table  9. 


TABLE  9.  MULTIPLE-COMPONENT  TEST  MATRIX 


Chemical 

Target  Concentrations. 

mq/m3  for 

Runs 

1 

2 

3 

4 

5 

6 

7 

8 

DCE 

0 

10 

0 

0 

0 

0 

0 

0 

TCE 

0 

50 

15 

0 

0 

0 

0 

0 

Benzene 

15 

0 

40 

800 

50 

300 

300 

30 

Toluene 

30 

0 

0 

100 

30 

300 

300 

300 

Et-benzene 

100 

0 

0 

0 

10 

300 

30 

30 

o-Xylene 

0 

0 

0 

0 

0 

300 

30 

30 

Results  for  these  runs  are  summarized  in  Table  10. 


These  data  lead  to  two  conclusions: 


1.  The  breakthrough  times  are  reduced  considerably  when  more  than  one 
component  is  present,  with  the  more  volatile  components  being 
affected  most  severely. 

2.  The  DRW  model  predictions,  while  not  perfect,  are  sufficiently 
accurate  to  provide  a  basis  for  estimating  behavior  of  other 
multiple-component  systems  (e.g.,  other  combinations  of  sorbates 
and/or  concentrations). 

Because  the  slower  moving  zones  in  the  carbon  bed  at  least  partially 
desorb  the  leading  zones,  the  gas-phase  concentrations  of  various 


TABLE  10.  MULTICOMPONENT  BREAKTHROUGH  TIMES 


(Sorbent  =  0.134  g  CT) 


Sorbate 

C,  ing/m^ 

t/t0* 

1  tp/t0b 

Run  No.  1 

:  F  = 

165  cc/min 

Benzene 

15 

0.18 

0.10 

Toluene 

50 

0.40 

0.35 

Ethylbenzene 

111 

0.94 

0.94 

Run  No.  2 

:  F  = 

166  cc/min 

1 , 1-Dichl oroethy 1 ene 

11 

0.62 

0.63 

Trichloroethylene 

49 

1.01 

1.00 

Run  No.  3 

:  F  = 

122  cc/min 

Trichloroethylene 

11 

0.20 

0.19 

Benzene 

37 

0.72 

0.88 

Run  No.  4 

:  F  = 

159  cc/min 

Benzene 

802 

0.81 

0.83 

Toluene 

88 

0.54 

0.56 

Run  No.  5 

:  F  * 

159  cc/min 

Benzene 

49 

0.57 

0.56 

Toluene 

29 

0.69 

0.68 

Ethylbenzene 

8 

0.53 

0.69 

Run  No.  6a 

:  F  = 

159  cc/min 

Benzene 

350 

0.38 

0.13 

Toluene 

344 

0.41 

0.21 

Ethylbenzene 

308 

0.47 

0.34 

o-Xylene 

304 

0.61 

0.56 

Run  No.  6b 

:  F  = 

165  cc/min 

Benzene 

372 

0.44 

0.18 

Toluene 

352 

0.53 

0.28 

Ethylbenzene 

321 

0.70 

0.48 

o-Xylene 

29 

0.17 

0.12 

23 


TABLE  10.  MULTICOMPONENT  BREAKTHROUGH  TIMES  (CONCLUDED) 
(Sorbent  =  0.134  g  CT) 


Sorbate 

C,  mg/m3 

t/t0a 

tp/t0b 

Benzene 

341 

0.48 

0.32 

Toluene 

309 

0.70 

0.50 

Ethylbenzene 

40 

0.32 

0.23 

o-Xylene 

37 

0.42 

0.39 

Run  No. 

8  :  F  = 

169  cc/min 

Benzene 

30 

0.14 

0.07 

Toluene 

332 

0.81 

0.58 

Ethylbenzene 

32 

0.27 

0.21 

o-Xylene 

30 

0.38 

0.38 

aRatio  of  observed  10  percent  breakthrough  time  to 
predicted  single- component  breakthrough  time. 


Predicted  ratio  of  multicomponent  breakthrough  time 
to  single-component  breakthrough  time. 

components  in  different  parts  of  the  bed  are  not  necessarily  the  same  as 
those  in  the  feedstock.  This  is  readily  obvious  from  observation  of  the 
eluent  concentrations  in  these  multiple-component  experiments.  Figure  5 
shows  the  eluent  concentrations  measured  for  Run  6  of  this  series  of 
experiments.  The  concentrations  of  each  of  the  three  faster-moving  compo¬ 
nents  go  through  maxima  that  appreciably  exceed  their  respective  input 
concentrations.  This  aspect  of  multicomponent  behavior  is  qualitatively 
predicted  by  the  DRW  model.  For  example,  Figure  6  shows  the  calculated 
concentrations  for  the  adsorbed  phase  in  each  bed  zone  for  Run  6  although 
these  are  shown  as  smooth  curves  drawn  through  the  calculated  points,  and 
are  probably  not  very  accurate,  the  suggested  trends  are  significant: 

1.  Benzene  is  present  to  some  extent  in  all  four  zones,  but  It  forms 
a  significant  fraction  of  the  adsorbed  phase  only  in  the  last 
zone. 


4  Component  Run  with  CT  as  the  Sorbent 


Distribution  of  Components  in  Adsorbed  Phase 


Figure  6.  Calculated  Adsorbed  Phase  Compositions  for  Run  Number 


2.  Toluene  dominates  zone  3;  ethylbenzene  dominates  zone  2,  but  is 
also  present  in  a  significant  portion  of  zone  1. 

Thus,  by  the  time  zone  4  reaches  the  bed  exit,  it  contains  almost  all 
of  the  benzene  adsorbed  during  the  run,  and  its  release  results  in  a  con¬ 
centration  in  the  gas  phase  that  greatly  exceeds  the  input  level.  In  this 
run  we  see  examples  that  suggest  both  extremes  noted  above: 

1.  Ethylbenzene  and  xylene  coadsorb  and  exclude  most  of  the  more 
volatile  components. 

2.  Benzene  is  largely  excluded  from  coadsorption  with  the  less 
volatile  components. 

Note: 

These  are  not  general  conclusions,  but  depend  on  the  particular  com¬ 
bination  of  sorbates  and  concentrations. 

E.  RELATIVE  HUMIDITY  EFFECTS 

Water  adsorption  isotherms  (Figures  B-l  to  B-8),  show  that  pore  con¬ 
densation  of  water  on  carbons  generally  does  not  occur  at  relative 
humidities  below  about  40-50  percent.  Therefore,  the  use  of  a  pseudo- 
Dubinin-Radushkevlch  Isotherm  to  represent  water  adsorption  offers,  at 
best,  a  crude  approximation.  Furthermore,  because  most  VOCs  are  only 
slightly  soluble  in  water,  the  use  of  ideal  solution  theory  with  such 
solutions  is  inappropriate.  Therefore,  one  cannot  use  the  DRW  model 
developed  on  this  program  to  estimate  relative  humidity  effects  with  more 
than  semi qualitative  accuracy.  It  has  been  suggested  (Werner)  that  for 
dual -component  sorbate  systems  in  which  one  of  the  components  is  water,  a 
pseudo-DR  Isotherm  can  be  used,  along  with  the  known  water  solubility  of 
the  second  component.  Such  an  approach  is  potentially  useful  only  if 
sufficient  Information  on  water  solubilities  is  available.  Unfortunately, 
water  solubilities  for  mixtures  of  VOCs  are  often  complex  functions  of  VOC 
concentrations,  and  complete  experimental  descriptions  of  systems  of 
Interest  are  not  generally  available.  Nonetheless,  it  is  useful  to 
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consider  the  qualitative  predictions  of  the  DRW  model  with  respect  to  the 
effect  of  relative  humidity  on  VOC  adsorption. 

If  we  estimate  the  affinity  coefficient  of  water  using  the  molar 
refractivlty  correlation  (see  Appendix  A),  we  find  that  the  affinity 
coefficient  for  water  is  0.24  times  that  of  benzene  or  trichloroethylene 
and  0.15  times  that  of  o-xylene  or  ethylbenzene.  Any  two  sorbates  compared 
will  have  equal  adsorption  potentials  when  the  ratios  of  the  logarithms  of 
their  relative  pressures  to  their  affinity  coefficients  are  equal  (assuming 
a  DR  isotherm).  We  can  use  this  to  estimate  the  relative  pressure  of  VOC 
needed  to  compete  on  an  equal  basis  with  water  vapor  as  function  of  the 
relative  pressure  of  water,  as  shown  in  Figure  7. 

The  curves  shown  in  Figure  7  (which  are  Independent  of  the  carbon) 
emphasize  the  fact  that  the  adsorption  potential  of  a  given  sorbate  Is  a 
function  of:  (a)  the  equilibrium  vapor  pressure  of  the  sorbate;  (b)  the 
actual  vapor  concentration;  and  (c)  the  affinity  coefficient  of  the 
sorbate.  At  high  values  of  P/P0,  almost  any  VOC  would  be  expected  to  be 
relatively  unaffected  by  the  presence  of  water  vapor.  At  low  P/P0  values, 
sorbates  having  low  affinity  coefficients  are  most  strongly  affected  by 
water  vapor.  The  curves  also  illustrate  the  magnitude  of  the  displacement 
effects  with  different  VOCs  at  low  concentrations.  For  example,  benzene  at 
P/P0  =  3  x  10~6  (about  1  mg/m3  at  25  C)  would  be  displaced  by  xylene  at 
P/P0  >  2  x  10-9  (about  6.4  x  10"5  mg/m3). 

General  conclusions  concerning  the  effect  of  water  vapor  that  can  be 
drawn  from  the  DRW  model  are: 

1.  The  higher  the  adsorption  potential  for  a  particular  component  the 
less  It  will  be  affected  by  the  presence  of  other  sorbates, 
including  water  vapor. 

2.  Sorbates  having  relatively  low  affinity  coefficients  should  show 
greater  concentration  effects  in  the  presence  of  other  sorbates. 
This  means  that  concentration  effects  with  o-xylene  and  water 
vapor  should  be  less  than  those  observed  by  Werner  for 
trichloroethylene  and  water  vapor. 

3.  There  should  be  no  differences  between  behaviors  of  different 
sorbates  when  comparison  is  made  on  the  basis  of  equal  adsorption 
potentials  rather  than  on  an  equal  concentration  basis. 
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-Bonzono  or  TCCi  y-Xylono  or  Ethylbonzone 


Relative  Humidity.  Par cant 

Figure  7.  Concentration  of  VOC  Required  for  Equal  Competition  with  Water  Vapor. 
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Initial  experiments  to  study  the  effect  of  water  vapor  on  VOC 
adsorption  on  carbons  consisted  of  replication  of  selected  experiments 
performed  by  Werner  on  an  earlier  program  for  the  Air  Force.  For  these 
'  experiments,  samples  of  the  carbon  used  by  Werner  were  obtained  and  were 

|  crushed  and  sieved  to  60/80  mesh.  A  comparison  of  Werner's  results  with 

|  trichloroethylene  adsorption  and  those  obtained  in  the  current  work  is 

given  In  Figure  8.  This  figure  shows  that  relatively  good  agreement  was 
obtained  with  Werner's  results. 

Additional  experiments  were  conducted  at  various  humidities  with 

1.1- dichloroethylene,  benzene,  and  o-xylene  as  the  sorbates  and  CT  as  the 
sorbent.  Data  from  these  experiments  are  shown  in  Figures  9  and  10. 

While  there  Is  somewhat  more  scatter  In  the  data  for  benzene  and 

1.1- dichloroethylene  than  for  o-xylene,  it  can  be  seen  that  1,1-dichloro¬ 
ethylene  and  benzene  are  more  severely  affected  by  the  presence  of  water 
vapor,  and  that  significant  effects  on  the  adsorption  of  1,1-dichloro¬ 
ethylene  and  benzene  occur  at  humidities  above  about  30  percent.  Further¬ 
more,  the  lower  the  concentration  of  1,1-dichloroethylene,  the  more  severe 
the  humidity  effect.  This  agrees  with  Werner's  findings  and  the  general 
Inferences  based  on  the  DRW  model.  Experiments  were  also  conducted  to 
determine  the  synergistic  effects  of  temperature  and  humidity  of  the 
adsorption  of  mixtures  of  VOCs.  In  these  experiments,  the  feedstock  was 
humidified  to  approximately  75  percent  RH  at  25  degrees  Centigrade,  while 
the  bed  was  operated  at  various  temperatures  from  25  to  55  degrees 
Centigrade.  The  net  effect  of  Increased  temperature  In  the  bed  Is  to  lower 
the  effective  relative  humidity  within  the  bed.  Results  of  these  runs  are 
given  in  Table  11.  In  this  table,  breakthrough  times  are  given  as  ratios 
of  measured  breakthrough  times  to  those  calculated  for  the  dry  system. 

This  normalization  procedure  removes  some  of  the  variation  In  the  data  due 
to  slightly  different  conditions  used  for  the  individual  experiments.  For 
convenience  In  Interpretation,  these  data  are  plotted  in  Figure  11,  which 
shows  that  xylene  Is  affected  least  by  the  high  humidity,  and  the  two 
chlorinated  compounds  are  affected  most  strongly.  This  is  in  agreement 
with  the  relative  adsorption  potentials  of  these  compounds.  All  five 
compounds  show  Increased  retention  with  Increased  temperatures  to  35-45 
degrees  Centigrade. 
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elative  Humidity  on  1,1-Dlchloroethylene  Adsorption 


TABLE  11.  INITIAL  HUMIDITY/TEMPERATURE  RUNS 
(CARBON  BED  =  CT) 


Sorbate 


C,  mg/m3  t10tio(cal) 


Run 

1111;  F=150 

cc/min; 

T 

= 

26  C;  RH 

=  743 

Benzene 

356 

0.67 

Toluene 

328 

0.79 

o-Xylene 

337 

0.85 

Run 

1112;  F= 1 52 

cc/mln; 

T 

= 

35  C;  RH 

=  75 

Benzene 

346 

1.22 

Toluene 

328 

1.01 

o-Xylene 

336 

0.87 

Run 

1113;  F=146 

cc/min; 

T 

35  C;  RH 

=  69 

Benzene 

350 

1.22 

Toluene 

332 

0.99 

o-Xylene 

345 

0.89 

Run 

1114;  F=143 

cc/mln; 

T 

3 

45  C;  RH 

=  73 

Benzene 

343 

1.09 

Toluene 

326 

0.81 

o-Xylene 

340 

0.82 

Run 

1125;  F=146 

cc/min; 

T 

= 

55  C;  RH 

=  75 

Benzene 

223 

109 

• 

Toluene 

205 

0.93 

o-Xylene 

156 

0.82 

Run 

1118;  F= 1 48 

cc/min; 

T 

= 

25  C;  RH 

=  70 

1,1-Dichloroethylene 

304 

0.17 

Trichloroethylene 

280 

0.35 

Run 

1119;  F= 1 50 

cc/min; 

T 

3 

25  C;  RH 

=  65 

1,1-Dichloroethylene 

258 

0.21 

Trichloroethylene 

258 

0.36 

Run 

1120;  F= 1 50 

cc/min; 

T 

= 

36  C;  RH 

=  69 

1,1-Dichloroethylene 

259 

Trichloroethylene 

262 

0.78 

Run 

1121;  F  =  1 4  3 

cc/min; 

T 

- 

45  C;  RH 

=  66 

1 , 1-Dichloroethy lene 

283 

0.84 

Trichloroethylene 

279 

0.75 

1 , 1  -Diehl oroethy iene 
Trichloroethylene 


TABLE  11.  INITIAL  HUMIDITY/TEMPERATURE  RUNS 
(CARBON  BED  =  CT)  (CONCLUDED) 


Sorbate 


C,  mg/m3  tiotio(cal) 


Run  1122;  F= 1 48  cc/min;  T  =  25  C;  RH  =  5 
1,1-Dichloroethylene  293  1.01 

Trichloroethylene  307  0.96 

Run  1124;  F= 150  cc/min;  T=  35  C;  RH  =  72 
317  0.81 

277  0.80 


1,1-Dichloroethylene 

Trichloroethylene 


Increased  Bed  Temperature  on  Retention  of  VOCs  from  Humidified  Gas  Streams 


SECTION  IV 
ECONOMIC  ASSESSMENT 


The  objective  of  this  subtask  was  to  provide  a  preliminary  economic 
assessment  of  various  carbon  adsorption  systems  for  treatment  of  well  water 
contaminated  with  volatile  organic  compounds  (VOCs).  The  following  treat¬ 
ment  options  were  included  in  the  analysis: 


Option  1A:  Air  stripping  of  well  water,  followed  by  air  cleaning  with 
activated  carbon  for  Case  1  contaminants  listed  in 
Table  12. 

Option  IB:  Air  stripping  of  well  water,  followed  by  air  cleaning  with 
activated  carbon  for  Case  2  contaminants  listed  in 
Table  12. 

Option  1C:  Air  stripping  of  well  water  followed  by  air  cleaning  with 
activated  carbon  for  Case  3  contaminants  listed  in 
Table  12. 

Option  2A:  Air  stripping  of  well  water,  followed  by  preheating  of  the 
contaminated  air  and  cleaning  of  the  heated  air  with 
activated  carbon  for  Case  1  contaminants. 

Option  2B:  Air  stripping  of  well  water,  followed  by  preheating  of  the 
contaminated  air  and  cleaning  of  the  heated  air  with 
activated  carbon  for  Case  2  contaminants. 

Option  2C:  Air  stripping  of  well  water  followed  by  preheating  of  the 
contaminated  air  and  subsequent  cleaning  of  the  heated  air 
with  activated  carbon  for  Case  3  contaminants. 

Option  3A:  Direct  contact  of  well  water  with  activated  carbon  for 
Case  1  contaminants. 

Option  3B:  Direct  contact  of  well  water  with  activated  carbon  for 
Case  2  contaminants. 

Option  3C:  Direct  contact  of  well  water  with  activated  carbon  for 
Case  3  contaminants. 


TABLE  12.  VOC  CONCENTRATION  IN  WELL  WATER 


VOC 

VOC 

Concentration, 

ppb 

Case  1 

Case  2 

Case  3 

Benzene 

100 

0 

1000 

Toluene 

100 

0 

1000 

Ethylbenzene 

100 

0 

1000 

Xylene 

100 

0 

1000 

Trichloroethylene 

0 

100 

0 

1, 1-Dichloroethylene 

0 

100 

0 

A.  PRELIMINARY  PROCESS  DESIGN 

A  preliminary  process  design  was  performed  to  size  major  process 
equipment  items  for  estimation  of  capital  costs.  Operating  costs  were 
estimated  from  utility  costs,  carbon  usage,  and  capital-related  costs,  such 
as  amortization  and  maintenance  costs.  A  process  design  was  based  on  the 
criteria  listed  below: 

a.  Well  water: 

Treatment  rate:  50  gpm  12  hours/day,  365  days/year 
Treatment  period:  10  years 

Temperature:  13  degrees  Centigrade  (55  degrees  Fahrenheit) 

b.  Stripping  air: 

Flow  rate:  335  cfm  (air-to-water  ratio  of  50) 

Temperature:  20  degrees  Centigrade  (68  degrees  Fahrenheit) 

Design  procedures  and  results  are  described  below. 
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B.  AIR-STRIPPING  TOWER 

Packed  towers  are  the  most  effective  and  widely  used  equipment  for  the 
aeration  process.  A  packed  tower  with  a  2-foot  diameter  and  22.5-foot 
packing  height  was  sized  on  discussions  with  an  equipment  supplier  (Refer¬ 
ence  6).  The  superficial  gas  velocity  in  this  tower  is  estimated  at  around 
1.8  ft/sec,  and  the  liquid  (water)  rate  at  around  8000  lb/hr-ft^.  The 
pressure  drop  through  1-inch  Raschig  ring  packing  is  estimated  at  around 
0.5  inch  HpO  per  foot  of  packing,  normally  found  at  the  flooding  point  and, 
therefore,  is  acceptable. 

C.  CARBON  ADSORPTION  TOWER 

An  adsorption  tower  for  air  cleaning  was  sized  at  4-foot  diameter  by 
4-foot  bed  depth,  which  holds  approximately  1500  pounds  of  activated 
carbon,  based  on  discussions  with  a  supplier  of  activated  carbon  adsorption 
systems  (Reference  7).  The  superficial  gas  velocity  for  this  unit,  based 
on  an  airflow  rate  of  335  cfm,  is  estimated  at  27  ft/min,  which  is  within 
the  range  of  20  to  100  ft/min  of  superficial  gas  velocity  utilized  normally 
in  vapor  phase  system  designs  (Reference  7).  The  bed  depth  of  4  feet  is 
conservative  when  compared  with  normal  bed  depths  in  the  2-  to  3-foot  range 
(Reference  7). 

For  a  liquid-phase  contact  system  to  be  utilized  for  direct  purifica¬ 
tion  of  the  contaminated  well  water  without  air  stripping,  a  disposable 
adsorption  system  utilizing  four  disposable  drum  adsorbers  also  was 
selected  (Reference  7).  The  drums  are  connected  in  two  parallel  lines  with 
two  drums  in  series  in  each  line.  Each  drum  contains  1000  pounds  of 
activated  carbon.  The  drum  Is  approximately  44  inches  in  diameter  and 
67  Inches  in  height.  The  drum  is  supplied  with  inlet  and  outlet  pipe 
connections  and  can  be  easily  installed  and  removed  from  the  process. 

After  the  carbon  in  the  drum  is  used  up,  the  drum  is  taken  off  the  process 
and  sent  to  a  hazardous  landfill  for  disposal. 
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0.  AIR  HEATER 


Assuming  that  the  air  stream  at  the  outlet  of  the  air  stripper  reaches 
the  inlet  water  temperature  of  55  degrees  Fahrenheit  (13  degrees  Centi¬ 
grade)  and  a  90  percent  relative  humidity,  the  partial  pressure  of  water 
vapor  in  the  air  stream  will  be  10.0  mm  Hg.  Published  data  by  Werner 
(Reference  8)  indicate  that  sorption  capacity  of  activated  carbon  levels 
off  at  around  25  percent  relative  humidity.  To  reduce  the  relative 
humidity  from  90  percent  to  25  percent,  the  air  stream  needs  to  be  heated 
from  55  degrees  Fahrenheit  (13  degrees  Centigrade)  to  93  degrees  Fahrenheit 
(34  degrees  Centigrade).  The  theoretical  energy  requirement  for  heating 
the  air  stream  is  estimated  at  around  15,000  Btu/hr,  or  4.1  kW  electric. 


E.  ADSORPTION  CAPACITY  OF  ACTIVATED  CARBON 


Based  on  a  complete  removal  of  VOCs  by  air  stripping,  the  V0C  concen¬ 
tration  in  the  stripping  air  is  estimated  at  1.65  ppm  by  weight  (2.2  mg/m3) 
for  each  of  the  Case  1  contaminants  (i.e.,  benzene,  toluene,  ethylbenzene, 
and  o-xylene),  and  for  each  of  the  Case  2  contaminants  (i.e.,  trichloro¬ 
ethylene  and  dichloroethylene) .  (Note:  This  concentration  is  approxi¬ 
mately  one  order  of  magnitude  lower  than  was  specified  for  the  experimental 
phase  of  this  work.)  For  Case  3  contaminants,  the  V0C  concentration  in  the 
stripping  air  is  estimated  at  16.5  ppm  by  weight  for  each  V0C.  Gas-phase 
adsorption  capacities  calculated  from  the  adsorption  test  results  obtained 
in  this  project  are  given  in  Table  13  for  the  Case  1,  Case  2,  and  Case  3 
contaminants  with  and  without  heating  of  the  stripper  air  (i.e.,  at 
55  degrees  Fahrenheit  and  90  percent  RH  without  heating  and  at  93  degrees 
Fahrenheit  and  25  percent  RH  with  heating).  Daily  carbon  usage  was  calcu¬ 
lated  from  a  daily  generation  rate  of  0.03  lb/day  (for  Case  1  and  Case  2 
contaminants)  and  0.3  lb/day  (for  Case  3  contaminants)  for  each  V0C  and  the 
adsorption  capacity  for  individual  V0C.  The  total  carbon  usage  for  a 
mixture  of  VOCs  was  obtained  as  the  sum  of  the  individual  carbon  require¬ 
ments.  The  useful  life  of  the  1500-pound  carbon  bed  for  vapor  phase  system 
was  calculated  from  the  total  daily  carbon  usage.  Calculated  results  are 
shown  in  Table  14.  Estimated  bed  lives  are  short  (i.e.,  less  than  10  days) 
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TABLE  13.  CALCULATED  GAS-PHASE  CARBON  ADSORPTION  CAPACITIES 


1 , 1 -0 1 ch 1 oroet hy 1 ene 

ppm  by  wt. 

Temp. , 

percent 

mg 

VOC/g  carbon 

VOC 

VOC  Cone., 
ppm  by  wt. 

Air 
Temp. , 

F 

RH, 

percent 

mg 

Adsorption 
Capacity, 
VOC/g  carbon 

Benzene 

1.65 

55 

90 

0.1 

Toluene 

1.65 

55 

90 

0.5 

Ethylbenzene 

1.65 

55 

90 

2.0 

o-Xylene 

1.65 

55 

90 

3.0 

Benzene 

1.65 

93 

25 

0.2 

Toluene 

1.65 

93 

25 

1.2 

Ethylbenzene 

1.65 

93 

25 

5.7 

o-Xylene 

1.65 

93 

25 

8.7 

Benzene 

16.5 

55 

90 

0.9 

Toluene 

16.5 

55 

90 

4.4 

Ethylbenzene 

16.5 

55 

90 

14.6 

o-Xylene 

16.5 

55 

90 

18.1 

Benzene 

16.5 

93 

25 

48.9 

Toluene 

16.5 

93 

25 

113.3 

Ethylbenzene 

16.5 

93 

25 

152.9 

o-Xylene 

16.5 

93 

25 

112.0 

1,1-Dichloroethylene 

1.65 

55 

90 

20.0 

Trichloroethylene 

1.65 

55 

90 

0.1 

1,1-Dichloroethylene 

1.65 

93 

25 

40a 

Trichloroethylene 

1.65 

93 

25 

0.2 

Extrapolated  from  trichloroethylene  data. 


except  for  Case  3  contaminants  with  heating  of  stripper  air  with  an 
estimated  bed  life  of  112  days.  A  regenerable  system  is  the  more  economi¬ 
cal  choice  when  the  carbon-bed  life  is  expected  to  be  less  than  3  months 
(Reference  7).  For  the  purpose  of  the  present  economic  analysis,  it  was 
assumed  that  the  carbon  bed  for  vapor-phase  adsorption  is  regenerated  in 
all  cases.  During  regeneration,  steam  is  passed  through  the  carbon  bed  to 
desorb  the  VOCs  and  collected  outside  the  bed  by  condensation.  The  steam 
requirement  Is  around  1.5  lb/lb  carbon  (i.e.,  around  2250  pounds  of 
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TABLE  14.  ESTIMATED  CARBON-BED  LIFE  FOR  GAS-PHASE  ADSORBER 


voc 

Daily  Carbon 
Usage,  lb/day 

Estimated  Beda 
Life,  day 

Without 

Heating 

With 

Heating 

Without 

Heating 

With 

Heating 

Case  1  Contaminants 

Benzene 

300 

150.0 

Toluene 

60 

25.0 

Ethylbenzene 

15 

5.3 

o-Xylene 

10 

3.4 

TOTAL 

385 

184 

3.9 

8.2 

Case  2  Contaminants 

1,1-Dlchloroethylene 

1.5 

0.75 

Trichloroethylene 

300 

150 

TOTAL 

302 

151 

5.0 

9.9 

Case  3  Contaminants 

Benzene 

333 

6.13 

Toluene 

68 

2.65 

Ethylbenzene 

21 

1.96 

o-Xylene 

17 

2.67 

TOTAL 

439 

13.4 

3.4 

112 

aBased  on  1500-pound  carbon  bed. 


steam  will  be  used  to  regenerate  the  1500-pound  carbon  bed).  The  steam 
condensate  will  be  contaminated  with  the  VOCs  and  must  be  disposed  of 
safely.  The  condensate  can  be  disposed  of  by  either:  (1)  the  liquid  can 
be  hauled  away  directly,  or  (2)  the  liquid  can  be  further  treated  by  carbon 
adsorption  to  reduce  the  waste  volume,  and  the  resulting  spent  carbon  can 
be  hauled  away. 

Liquid-phase  adsorption  capacities  were  obtained  from  a  manufacturer 
of  activated  carbons.  Adsorption  capacity  data  are  given  in  Table  15  for 
Case  1  and  Case  2  contaminants  at  0.1  mg/L  (100  ppb)  for  direct  treatment 
of  contaminated  well  waters  and  at  10  mg/L  (10  ppm)  for  treatment  of  steam 
condensate  resulting  from  regeneration  of  the  vapor-phase  adsorption 


TABLE  15.  LIQUID-PHASE  CARBON  ADSORPTION  CAPACITIES 
AT  25  DEGREES  CENTIGRADE  (Reference  9) 


Liquid-Phase 

Concentration, 

mg/L 


Adsorption 
Capacity, 
mg  VOC/g  carbon 


a0ata  for  p-xylene. 


system,  and  Case  3  contaminants  at  1  mg/L  (1000  ppb)  for  direct  treatment 
of  contaminated  well  waters  and  at  100  mg/L  (100  ppm)  for  treatment  of 
steam  condensate  resulting  from  regeneration  of  the  vapor-phase  adsorption 
system. 

Daily  carbon  usage  was  calculated  from  a  daily  generation  rate 
of  0.03  lb/day  for  each  VOC  of  Case  1  and  Case  2  contaminants  and  a  daily 
generation  rate  of  0.3  lb/day  for  each  VOC  of  Case  3  contaminants,  and  the 
individual  adsorption  capacity.  The  total  carbon  usage  for  a  VOC  mixture 


was  obtained  as  the  sum  of  the  carbon  usages  for  the  individual  VOCs.  The 
useful  life  of  the  disposable  drum  adsorber,  containing  1000  pounds  of  car¬ 
bon,  was  calculated  from  the  total  daily  carbon  usage.  Results  shown  in 
Table  16  indicate  that  a  1000-pound  carbon  bed  would  last  from  9  to  498 
days,  depending  on  the  VOC  types,  and  the  VOC  concentrations  in  the  liquid 
phase  (i.e.,  well  water  or  steam  condensate). 

TABLE  16.  ESTIMATED  CARBON-BED  LIFE  FOR  LIQUID-PHASE  ADSORBER 


Daily  Carbon 

Estimated  Bedb 

Usaqe, 

Ib/day 

Life, 

day 

0.1 

10 

0.1 

10 

VOC 

mg/La 

mg/La 

mg/La 

mg/La 

Case  1  Contaminants 
Benzene 

21.4 

4.17 

Toluene 

3.3 

0.42 

Ethylbenzene 

3.3 

0.10 

o-Xylene 

0.5 

0.25 

TOTAL 

28.5 

4.94 

38.8 

202 

Case  2  Contaminants 
1,1-01  chi oroethylene 

21.4 

1.76 

Trichloroethylene 

4.4 

0.25 

TOTAL 

25.8 

2.01 

38.8 

498 

Case  3  Contaminants 
Benzene 

88.2 

3.0 

Toluene 

12.0 

1.5 

Ethylbenzene 

5.7 

0.4 

o-Xylene 

3.6 

1.5 

TOTAL 

109.5 

33.4 

9.1 

29.9 

aV0C  concentration  in  the  liquid  phase. 
bBased  on  1000-pound  carbon  bed. 


The  processes  selected  and  described  above  for  the  various  treatment 
options  can  be  summarized  as  follows: 


Options  1A,  IB,  and  1C:  The  well  water  is  decontaminated  using  a 
packed-bed,  air-stripping  tower.  The  off-gas  from  the  stripping  tower 
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containing  VOCs,  is  purified  using  a  regenerable  granular  carbon 
adsorber.  After  the  carbon  in  the  adsorber  is  saturated,  it  is 
regenerated  with  steam.  The  regeneration  requires  approximately 
1  hour,  and  can  be  scheduled  during  the  12-hour  period  when  the  air 
stripper  is  down.  The  steam  condensate  collected  from  carbon  regen¬ 
eration  is  treated  with  a  disposable  drum  carbon  adsorber.  The  con¬ 
densate  treatment  system  consists  of  two  disposable  units  in  series. 
After  usage,  the  disposable  drum  containing  saturated  carbon  is 
replaced  with  a  fresh  unit,  and  the  used  unit  is  hauled  away  for 
offsite  disposal  at  a  chemical  waste  landfill. 

Options  2A,  2B,  and  2C:  The  process  steps  for  Options  2A,  2B,  and  2C 
are  the  same  as  for  Options  1A,  IB,  and  1C,  except  for  an  additional 
step  in  which  the  off-gas  stream  from  the  air  stripper  is  preheated 
before  it  enters  the  carbon  adsorber. 

Options  3A,  3B,  and  3C:  The  well  water  is  decontaminated  by  direct 
contact  with  granular  activated  carbon.  The  liquid-phase  adsorption 
system  consists  of  four  disposable  drum  carbon  adsorber  units,  con¬ 
nected  in  two  parallel  lines  with  two  units  in  series  in  each  line. 
After  usage,  the  disposable  drum  containing  saturated  carbon  is 
replaced  with  a  fresh  unit,  and  the  used  unit  is  hauled  away  for 
disposal  at  a  chemical  waste  landfill. 

F.  PROCESS  ECONOMICS 

Capital  and  operating  costs  were  derived  from  various  sources, 
Including  budgetary  estimates  provided  by  a  manufacturer  of  activated 
carbon,  equipment  vendors,  and  published  data.  Operating  costs  were 
estimated  with  an  assumption  that  the  facility  will  be  operated  unattended, 
except  for  one  operator  onsite  during  regeneration  of  carbon  adsorption 
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system  and  during  replacement  of  disposable  carbon  adsorbers.  The  operat¬ 
ing  costs  Included  the  following  Items: 

1.  Amortized  capital  cost  at  10  percent  of  capital  cost  per  year. 

2.  Maintenance  cost,  including  materials  and  labor,  at  7  percent  of 
capital  cost  per  year. 

3.  Utilities  (electricity,  steam). 

4.  Cost  of  activated  carbon  replacement  for  non-regenerable  adsorp¬ 
tion  systems. 

5.  Operator  during  adsorber  regeneration  and  replacement. 

Estimated  capital  and  operating  costs  are  summarized  In  Tables  17  and 
18,  respectively.  Capital  costs  are  identical  for  Case  1,  Case  2,  and 
Case  3  VOCs.  Capital  costs  for  the  liquid-phase  adsorption  system  without 
air  stripping  ( 1 . e . ,  Options  3A,  3B,  and  3C)  are  about  one-fourth  of  the 
capital  costs  for  the  vapor-phase  adsorption  systems  that  include  air 
stripping.  The  lower  capital  costs  are  due  to  a  relatively  simple  design 
based  on  the  use  of  disposable  carbon  adsorber  units. 

Operating  costs  for  the  liquid-phase  adsorption  system  are  about  20  to 
130  percent  higher  than  the  costs  for  the  vapor-phase  adsorption  system. 
Therefore,  we  would  recommend  the  vapor-phase  adsorption  system  ( 1 . e . ,  air 
stripping  followed  by  purification  of  the  stripper  air  by  carbon  adsorp¬ 
tion)  as  the  more  economical  choice.  Comparison  of  operating  costs  for 
Options  1A  and  IB  without  heating  of  the  stripper  air  and  for  Options  2A 
and  2B  with  stripper  air  heating  shows  about  10  percent  reduction  in  cost 
by  heating  the  stripper  air  for  Case  1  and  Case  2  contaminants.  However, 
there  Is  considerable  uncertainty  in  this  estimate  because  of  uncertainties 
in  the  adsorption  capacities  at  the  low  concentration  levels  that  were 
considered  in  the  economic  analysis.  Therefore,  we  cannot  conclusively 
state  that  this  estimated  economic  benefit  of  stripper  air  heating  Is 
significant  for  either  Case  1  or  Case  2  contaminants.  Comparison  of 
operating  costs  for  Option  1C  without  heating  of  the  stripper  air  and  for 
Option  2C  with  stripper  air  heating  shows  about  23  percent  reduction  In 
cost  by  heating  the  stripper  air  for  Case  3  contaminants.  Heating  of  the 
stripper  air  is,  therefore,  recommended  for  Case  3  contaminants. 


TABLE  17.  ESTIMATED  CAPITAL  COSTS 


1A,  IB,  1C  2A,  2B,  2C  3A,  3B,  3C 


Equipment  Cost,  $ 

Air-stripping  12,500  12,500 

Tower  (Reference  6) 

Water  Pump  (Reference  10)  500  500  500 

Air  Heater  (Reference  11)  --  1,000 

Carbon  Adsorber  36,000  36,000 

(air)  (Reference  7) 

Carbon  Adsorb er  6,400  6,400 

(condensate)  (Reference  9) 

Carbon  Adsorb er  --  --  12,800 

(well  water)  (Reference  9) 

Total  Purchased  55,400  56,400  13,300 

Equipment  Cost,  $ 

Capital  Costb,  $  105,000  107,000  25,000 

treatment  Options: 


1A  =  Air  stripping  of  water  followed  by  carbon  treatment  of 
air  (Case  1). 

IB  =  Air  stripping  of  water  followed  by  carbon  treatment  of 
air  (Case  2). 

1C  =  Air  stripping  of  water  followed  by  carbon  treatment  of 
air  (Case  3). 

2A  =  Air  stripping  of  water  followed  by  preheating  of  air  feed 
to  carbon  bed  (Case  1). 

2B  =  Air  stripping  of  water  followed  by  preheating  of  air  feed 
to  carbon  bed  (Case  2). 

2C  =  Air  stripping  of  water  followed  by  preheating  of  air  feed 
to  carbon  bed  (Case  3). 

3A  =  Direct  contact  of  water  with  carbon  (Case  1). 

3B  =  Direct  contact  of  water  with  carbon  (Case  2). 

3C  =  Direct  contact  of  water  with  carbon  (Case  3). 

^Capital  cost  was  obtained  by  applying  a  multiplication  factor 

of  1.89  to  the  total  purchased  equipment  cost  (Reference  12). 
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ESTIMATED  OPERATING  COSTS 


TABLE  18.  ESTIMATED  OPERATING  COSTS  (CONCLUDED) 


frVASftaeSe 


The  main  uncertainty  In  the  cost  estimates  Is  associated  with  extrap¬ 
olation  of  vapor-phase  adsorption  data  that  have  been  experimentally 
determined  at  relatively  high  VOC  concentrations  to  extremely  low  concen¬ 
trations  for  the  purpose  of  process  design  and  economic  analysis.  Addi¬ 
tional  vapor-phase  adsorption  data  are  needed  at  VOC  concentrations  below 
10  ppm  for  multicomponent  systems  that  would  be  directly  applicable  to 
process  design.  The  data  would  be  needed  to  refine  process  economics  and 
enable  the  selection  of  an  optimum  treatment  process  for  each  application. 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


Of  the  various  carbon  sorbents  considered  on  this  program,  the 
specialty  carbons  showed  no  special  advantages  over  the  more  standard 
activated  carbons,  and  would  prove  much  less  desirable  than  the  standard 
carbons  because  of  cost  factors  involved.  Within  the  group  of  standard 
carbons,  the  CT  and  SK-4  were  superior  to  LCL  and  WV-B  with  respect  to 
their  adsorption  capacities.  CT  and  SK-4  are  both  coconut-derived  carbons 
produced  by  the  same  manufacturer.  Differences  between  these  two  carbons 
may  be  because  they  are  produced  from  different  parts  of  the  coconut. 

The  model  based  on  the  Dubinin-Radushkevich  isotherm  and  the  Wheeler 
bed  permeation  equation  has  proven  valuable  for  correlating  adsorption  of 
VOCs,  and  is  especially  useful  for  prediction  of  multiple  component  adsorp¬ 
tion  behavior.  Even  though  the  DR  isotherm  is  not  strictly  applicable  to 
the  adsorption  of  water  vapor,  the  model  is  an  aid  to  understanding  the 
qualitative  aspects  of  the  effects  of  water  vapor  on  VOC  adsorption.  It  is 
recommended  that:  (1)  the  model,  in  its  current  form,  be  utilized  for 
estimation  of  bed  requirements  for  adsorption  of  complex  VOC  mixtures,  and 
(2)  that  consideration  be  given  to  refinement  of  the  model  with  respect  to 
relative  humidity  effects.  Such  refinement  will  require  incorporation  of  a 
data  base  of  solubilities  for  both  VOCs  in  water  and  water  in  VOC  mixtures, 
and  improved  simulation  of  water  isotherms. 

Results  of  the  current  work  indicate  that  significant  improvement  in 
retention  of  VOCs,  especially  for  the  more  volatile  VOCs,  can  be  obtained 
by  reducing  the  effective  relative  humidity  to  below  40-50  percent.  One 
means  for  accomplishing  this  is  by  increasing  the  bed  temperature.  Such 
increased  temperature  has  the  effect  of  decreasing  the  "dry  capacities"  of 
the  VOCs  by  virtue  of  the  fact  that  the  P0  values  are  increased.  However, 
the  temperature  coefficient  for  the  P0  of  water  is  greater  than  those  of 
the  VOCs  examined  on  this  program.  Thus,  increases  in  bed  temperature 
reduce  the  water  capacity  faster  than  they  reduce  the  VOC  capacities.  The 
net  result  is,  therefore,  a  gain  in  VOC  capacity  relative  to  the  "wet" 
system. 
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A  preliminary  process  design  and  economic  analysis  based  on  the  data 
obtained  from  the  experimental  program  indicates  that  air  stripping  of  the 
contaminated  well  water  followed  by  purification  of  the  stripper  air  by 
carbon  adsorption  is  more  economical  than  direct  contact  of  the  well  water 
with  a  carbon  adsorber.  An  economic  benefit  of  heating  the  stripper  air  is 
suggested  for  Case  3  contaminants  but  less  certain  for  either  Case  1  or 
Case  3  contaminants.  A  more  refined  cost  analysis  based  on  additional 
vapor-phase  adsorption  data  would  be  needed  to  further  optimize  the  treat- 
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DRW  MODEL 


for 

MULTICOMPONENT  ADSORPTION  OF  VOLATILE 
ORGANIC  HYDROCARBONS  ON  CARBON  SORBENTS 

A.  WORKING  MODEL  FOR  MULTICOMPONENT  BED  PERMEATION 

Considerable  work  has  been  done  in  recent  years  to  develop  methods  for 
correlation  and  prediction  of  adsorption  and  permeation  of  volatile  organic 
hydrocarbons  (VOC)  in  activated  carbon  beds.  An  excellent  review  of  this 
work  has  been  prepared  by  Werner  (M.  D.  Werner,  "Predicting  Gaseous  Phase 
Adsorption  by  Microporous  Adsorbents,"  CRC  Critical  Reviews  in  Environ¬ 
mental  Control,  16,  327-356,  1986).  Most  of  these  efforts  have  focused  on 
the  permeation  behavior  of  single  VOCs,  a  condition  that  is  rarely  repre¬ 
sentative  of  practical  sorbent  applications. 

With  the  more  practical  case  of  multiple  sorbates,  there  is  competi¬ 
tion  between  sorbates  in  the  adsorption  process  with  the  result  that  the 
overall  process  becomes  much  more  complex  to  treat  in  a  mathematically 
rigorous  fashion.  Nonetheless,  some  guidelines  are  needed  for  relating  the 
easily  accessible  single  component  data  to  multiple  component  behavior. 

The  following  is,  therefore,  offered  as  working  approximation  to  multicom- 
Donent  bed  permeation. 

B.  PHYSICAL  REPRESENTATION  AND  GLOSSARY  OF  TERMS 

For  purposes  of  this  discussion.  Figure  A-l  represents  a  sorbent  bed 
with  two  components,  a  and  b,  being  fed  to  the  bed  at  flow  rate  Q  and 
concentrations  Cao  and  Cbo-  Table  A-l  identifies  all  abbreviations  and 
symbols  to  be  used  in  this  discussion. 

Component  b  is  assumed  to  be  more  strongly  adsorbed  than  Component  a 
so  that  two  zones  of  differing  composition  develop  within  the  bed.  In 
zone  1,  both  components  are  present  at  gas  phase  concentrations  C^,  and 
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Figure  A-l.  Simplified  Multicomponent  Sorption  System. 

mole  fractions  Xu  in  the  adsorbed  phase.  Only  Component  a  is  present  in 
zone  2.  The  equilibrium  between  the  gas  phase  and  the  adsorbed  phase  for 
each  component  is  given  by  the  DR  isotherm,  viz. 


ei 

n—  =  exp  -  B 
woi 


RT  .  Poii 


in 


where  ideal  solution  behavior  is  assumed  for  the  adsorbed  phase  isee 
R.  J.  Grant  and  M.  Manes,  "Adsorption  Behavior  of  Binary  Hydrocarbon  Gas 
Mixtures  on  Activated  Carbon,"  I&EC  Fund.,  5,  490-498  (1966)1.  The 
breakthrough  time  for  zone  1  is  calculated  using  the  Wheeler  equation: 
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TABLE  A-l.  GLOSSARY  OF  SYMBOLS  AND  ABBREVIATIONS 


sorbate  index  jzone  index 

Adsorption  capacity 

g/g 

Saturation  capacity 

g/g. 

Gas  phase  concentration 

g/m3 

Bed  weight 

g. 

Gas  phase  volume 

m3 

Solid  phase  volume 

m3 

Flow  rate 

m3/min 

Sorbent  constant 

Temperature 

degrees  Kelvin 

Gas  constant  (1.987) 

cal/mole-deg 

Dubinin-Radushkevich 

DR  affinity  coefficient 

cal/mole 

Wheeler  rate  constant 

min-* 

Liquid  density 

g/m3 

Bulk  density  of  bed 

g/m3 

Density  of  sorbent 

g/m3 

Breakthrough  time  (10%) 

min 

Liquid  vapor  pressure 

atm 

Pressure  in  gas  phase 

atm 

Mole  fraction  in  adsorbed  phase 

Volume  fraction  in  adsorbed  phase 

where  the  component  having  the  longest  breakthrough  time  is  used  to  dis¬ 
tinguish  between  zones  1  and  2.  To  use  these  equations,  one  must  first 
determine  X^j.  To  do  this.  Grant  and  Manes  suggest  that  the  chemical 
potentials  for  all  components  in  a  given  zone  must  be  equal.  Therefore, 


Equations  3  and  4  can  generally  be  solved  by  the  Newton-Raphson  method 
to  yield  the  equilibrium  composition  of  the  adsorbed  phase.  To  proceed 
further,  it  is  necessary  to  know  the  values  of  B  and  Wo-jj.  Values  of  B  and 
Wo  are  usually  determined  using  a  reference  sorbate,  and  the  value  of  B  is 
constant  for  a  given  sorbent.  The  value  of  Wo  must  be  adjusted  for 
different  sorbates  because  of  differing  liquid  densities.  In  the  case  of 
multiple  sorbates,  care  must  be  taken  to  ensure  that  the  sum  of  all  of  the 
volumes  occupied  by  the  individual  sorbates  cannot  exceed  the  volume 
implied  by  the  value  of  Wo  for  the  reference  sorbate.  That  is. 


ref 
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where  the  values  of  Y-j j  are  derived  from  the  set  of  X^j  calculated  above. 
With  these  values  in  hand.  Equations  1  and  2  can  be  used  to  determine  the 
breakthrough  time  for  zone  1. 

The  computational  approach  for  the  successive  zones  is  similar,  but 
with  one  added  feature.  Solutions  for  multiple  zones,  and,  indeed,  experi¬ 
mental  evidence,  indicate  that  gas  phase  concentrations  in  the  successive 
zones  are  not  necessarily  the  same  as  the  concentrations  in  the  input 
stream.  If,  for  example.  Wen  is  ver-X  sma  1 1  and  Wejp  is  large,  then  the 
advancement  of  zone  1  will  necessarily  be  accompanied  by  desorption  of 
component  1.  This  will  have  the  effect  of  causing  Cn  to  be  greater  than 
Ci 2‘  Two  approaches  to  accommodate  this  effect  are: 

1.  C12  can  be  calculated  from  a  mass  balance  across  the  boundary 
between  zones  1  and  2. 

2.  The  advance  of  zone  1  can  be  viewed  in  terms  of  a  physical 
displacement  of  zone  2.  (This  is  equivalent  to  assuming  that  the 
isotherm  is  linear  between  Cn  and  C12O 

The  first  approach  is  preferable  because  of  the  fact  that  any  nonlinearity 
in  the  adsorption  isotherm  can  be  accommodated  through  an  iterative  compu¬ 
tation  of  the  bed  composition.  This  however  is  cumbersome,  and  the  solu¬ 
tion  converges  only  slowly.  On  the  other  hand,  the  second  approach  is 
direct.  It  is  equivalent  to  assuming  that  the  rates  of  movement  of  succes¬ 
sive  zones  are  additive,  i.e.. 
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The  overall  process  can,  thus,  be  summarized  in  terms  of  the  following 
computational  steps: 

1.  Calculate  the  composition  of  the  adsorbed  phase  in  the  first  zone 
using  Equations  3  and  4. 

2.  Use  Equations  1  and  2  to  determine  the  breakthrough  time  of  the 
slowest  moving  component. 

3.  Eliminate  the  slowest  moving  component  from  the  calculations  and 
repeat  step  1. 

4.  Use  Equations  1  and  6  to  calculate  the  breakthrough  time  for  the 
slowest  moving  component  of  zone  2. 

5.  Repeat  steps  3-5  for  successive  zones. 

Table  A-2  summarizes  the  results  of  application  of  this  simple  model 
in  comparison  with  experimentally  determined  breakthrough  times  observed 
with  selected  VOCs  on  an  activated  carbon.  In  this  case,  1 , 1-dichloroethy- 
lene  served  as  the  reference  sorbate  for  determining  the  values  of  B  and  Wo 
for  the  sorbent,  and  the  sorbent  was  an  activated  carbon  identified  as 
"CT."  In  general,  the  calculated  breakthrough  times  agree  reasonably  well 
with  the  observed  times  for  the  2  and  3  component  systems.  With  the 
4  component  systems,  the  calculated  results  are  qualitatively  correct,  but 
the  discrepancies  between  observed  and  calculated  times  are  greater, 
especially  for  the  more  rapidly  eluting  components.  Figure  A-2  shows  a 
typical  example  (Run  6a)  of  the  eluent  profiles  determined  by  periodic  gas 
chromatographic  analysis  at  the  bed  exit.  This  figure  clearly  illustrates 
the  separation  of  components  into  more  or  less  well-defined  zones. 
Furthermore,  concentrations  of  the  early  eluting  components  peak  out  at 
levels  considerably  higher  than  the  input  concentrations.  This  clearly 
indicates  the  effect  of  desorption  of  less  strongly  held  components. 
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TABLE  A-2.  MULTICOMPONENT  BREAKTHROUGH  TIMES 
(Sorbent  =  0.134  g  CT) 


Sorbate 

C,  mg/m^ 

t/toa  tp/tob 

Run  No.  1  : 

Q  « 

165  cc/min 

Benzene 

15 

0.18 

0.10 

Toluene 

50 

0.40 

0.35 

Ethylbenzene 

111 

0.94 

0.94 

Run  No.  2  : 

Q  = 

166  cc/min 

1,1-Dichloroethylene 

11 

0.62 

0.63 

Trichloroethylene 

49 

1.01 

1.00 

Run  No.  3  : 

Q  = 

122  cc/min 

Trichloroethylene 

11 

0.20 

0.19 

Benzene 

37 

0.72 

0.88 

Run  No.  4  : 

Q  - 

159  cc/min 

Benzene 

802 

0.81 

0.83 

Toluene 

88 

0.54 

0.56 

Run  No.  5  : 

Q  =  : 

159  cc/min 

Benzene 

49 

0.57 

0.56 

Toluene 

29 

0.69 

0.68 

Ethylbenzene 

8 

0.53 

0.69 

Run  No.  6a 

:  Q  = 

159  cc/min 

Benzene 

350 

0.38 

0.13 

Toluene 

344 

0.41 

0.21 

Ethylbenzene 

308 

0.47 

0.34 

o-Xylene 

304 

0.61 

0.56 

Run  No.  6b 

:  Q  - 

165  cc/min 

Benzene 

372 

0.44 

0.18 

Toluene 

352 

0.53 

0.28 

Ethylbenzene 

321 

0.70 

0.48 

o-Xylene 

29 

0.17 

0.12 

TABLE  A-2.  MULTICOMPONENT  BREAKTHROUGH  TIMES  (CONCLUDED) 
(Sorbent  =  0.134  g  CT) 


Sorbate 

C,  mg/m3 

t/toa 

tp/tob 

Run  No.  7 

:  Q  =  169 

cc/min 

Benzene 

341 

0.48 

0.32 

Toluene 

309 

0.70 

0.50 

Ethylbenzene 

40 

0.32 

0.23 

o-Xylene 

37 

0.42 

0.39 

Run  No.  8 

:  Q  =  169 

cc/min 

Benzene 

30 

0.14 

0.07 

Toluene 

332 

0.81 

0.58 

Ethylbenzene 

32 

0.27 

0.21 

o-Xylene 

30 

0.38 

0.38 

aRatio  of  observed  10  percent  breakthrough  time  to 
predicted  single  component  breakthrough  time. 


bRatio  of  predicted  multicomponent  breakthrough  time 
to  predicted  single  component  breakthrough  time. 


C.  ELUENT  CONCENTRATIONS 


This  simple  model  correctly  predicts  several  qualitative  aspects  of 
multiple  component  bed  permeation: 


1.  While  there  is  some  tendency  for  sharing  of  the  adsorbent  by 
multiple  components,  each  component  moves  through  the  bed  at 
different  rates.  This  results  in  the  separation  of  sorbates  and 
the  formation  of  zones  within  which  one  component  is  dominant. 

2.  The  advancement  of  the  slower  zones  tends  to  displace  the  zones 
formed  from  the  more  rapidly  moving  sorbates,  thus,  accelerating 
their  progress  through  the  bed. 

3.  The  ultimate  result  of  this  displacement  process  is  the 
consecutive  elimination  of  the  faster  moving  zones  from  the  bed. 
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Figure  A-2.  Eluent  Concentration  Profiles  for  Run  Number 


The  eluent  concentrations  can  be  estimated  as  follows: 

Ce  -j=C-jj+AC  (7) 

where  aC  is  obtained  from  a  mass  balance  at  the  bed  exit  -- 

QaC  =  Wei  i  (w/t i  i)(t-j_ii_i~ti  i)/tn)  (8) 

where  the  latter  factor  is  a  correction  for  the  fraction  of  the  bed 
actually  being  occupied  by  component  i  at  the  time  of  its  breakthrough. 

D.  RELATIVE  HUMIDITY  EFFECTS 

The  equations  given  above  assume  that  an  ideal  solution  is  formed 
between  the  sorbates  when  they  condense  in  the  sorbent  pores.  This  is  a 
reasonable  assumption  for  miscible  sorbates,  but  is  not  a  good  approxima¬ 
tion  for  cases  where  water  vapor  is  one  of  the  components.  Two  separate 
cases  should  be  considered  when  water  vapor  is  present: 

1.  At  high  relative  humidities,  water  will  tend  to  condense  in  the 
pore  structure  of  the  sorbent  and  the  saturation  pressure  of  the 
organic  vapor  is  PoX-j  where  Xj  is  the  liquid-liquid  solubility 
(expressed  as  a  mole  fraction)  of  the  organic  compound  in  water. 

2.  At  low  relative  humidities,  the  adsorbed  phase  is  the  organic  and 
the  saturation  pressure  of  the  water  vapor  is  given  by  PoXi  where 
X]  is  the  liquid-liquid  solubility  of  water  in  the  organic. 

In  neither  of  these  cases  can  the  solubilities  be  estimated  accurately 
by  ideal  solution  theory  or  any  other  simple  and  straightforward  process. 
Furthermore,  although  solubilities  of  many  individual  organic  compounds  in 
water  are  known,  data  on  the  solubilities  of  multiple  component  organic 
systems  in  water  and  data  on  the  solubility  of  water  in  organic  systems  are 
quite  limited.  In  spite  of  these  deficiencies,  the  DRW  model  still  has 
some  utility  for  estimating  the  effects  of  water  vapor  on  VOC  adsorption. 

The  approach  used  for  this  part  of  the  code  involves  simulation  of  the 
water  adsorption  isotherm  using  two  different  DR  isotherms:  one  set  of  DR 


parameters  is  used  for  the  low  humidity  side  of  the  Isotherm,  and  a  second 
set  of  parameters  is  used  at  the  high  humidities.  Figure  A-3  illustrates 
how  this  approach  approximates  the  water  adsorption  isotherm  of  CT.  The 
calculated  curve  does  not  have  exactly  the  same  shape  as  the  experimental 
curve.  Because  of  the  fact  that  the  intermediate  range  of  humidities 
(40~70  percent  RH)  are  of  least  interest  to  this  program,  the  curve  fitting 
was  done  so  as  to  best  simulate  the  experimental  data  below  40  percent  RH 
and  above  80  percent  RH.  This  approach  is  only  semiqual itatlve  at  best, 
but  even  qualitative  estimates  may  prove  useful  for  the  purposes  of  this 
program.  An  example  of  the  utility  of  this  approach  is  shown  in  Fig¬ 
ure  A-4.  This  figure  reproduces  the  data  shown  in  Figure  8,  with  an 
additional  curve  showing  the  calculated  effect  of  humidity  on  trichloro¬ 
ethylene  adsorption  on  CT  (DRW  parameters  for  Werner's  carbon  were  not 
available).  The  calculated  curve  shows  the  same  qualitative  behavior 
indicated  by  the  experimental  data,  i.e.,  relatively  little  effect  of  water 
vapor  at  humidities  less  than  30-40  percent,  and  a  strong  effect  at  higher 
humidities. 

E.  COMPUTER  CODE 

The  code  used  for  the  DRW  model  calculations  on  this  program  was 
written  in  Microsoft  Basic  for  use  on  IBM  compatible  systems.  Inasmuch  as 
the  development  of  this  code  was  not  one  of  the  principal  objectives  of  the 
current  work,  the  code  is  presented  merely  as  a  working  version  that  has 
not  been  optimized.  A  listing  of  the  code  is  given  following  this  text. 

All  necessary  sorbent  and  sorbate  parameters  are  included  In  the  code,  and 
the  only  input  required  is  the  specification  of  concentrations,  tempera¬ 
tures,  bed  weight,  and  flow  rate.  Units  required  are  specified  in  the 
prompts. 

10  REM  1 PREDICITI0N  OF  BREAKTHROUGH  CURVES  USING  WHEELER  AND  DUBININ 

RADUSHKEVICH  MODELS 

20  CIS: DIM  CI(10,10),X(10,10),T(10,10) 

30  OEFDBL  C,D,Y 

40  FOR  1*1  TO  7 .’READ  A(I),BP(I):NEXT 
50  FOR  1*1  TO  5: READ  W0( I ) ,K ( I ) : NEXT 
60  FOR  1-1  TO  5: READ  D(I) ,KR(I):NEXT 


Figure  A-3.  Calculated  Water  Adsorption  Isotherm 


70  FOR  1=1  TO  7:READ  B(I):NEXT 
80  FOR  1=1  TO  7 : READ  M(I):NEXT 
90  FOR  1=1  TO  7:READ  DL( I ) :NEXT 
100  FOR  1=1  TO  7: READ  NC$(I):NEXT 
110  FOR  1=1  TO  5:REA0  NS$(IP):NEXT 
120  PRINT  "Output  to  printer  ?  (Y/N)" 

130  A$=INKEY$:  IF  (A$="Y")  OR  (A$="y")  THEN  NFL=1  ELSE  IF  (A$="n")  OR 

(A$="N" )  THEN  NFL=0  ELSE  130 

140  FOR  1  =  1  TO  5 : PR  I  NT  I ;NS$( I ) : NEXT 

150  INPUT  "Select  one  of  the  above  sorbents  by  number  -  ",FS 
160  INPUT  "Number  of  sorbates  =  ",N 

170  CLS:PRINT  "Choose  Sorbates  from  the  following  list  by  responding  (Y/N)" 
180  PRINT:FOR  1=1  to  7:PRINT  NC$(I) 

190  A$= INKEYS: IF  A$="Y"  OR  A$="y"  THEN  FL(I)=1  ELSE  IF  A$="n"  OR 

A$="N"  THEN  FL(I)=0  ELSE  190 

200  IF  FL { I ) =1  THEN  K=K+1 

210  IF  K=N  THEN  220  ELSE  NEXT 

220  FOR  1=1  TO  6: IF  FL(I)  THEN  230  ELSE  240 

230  PRINT  "Concentration  (mg/m'3)  of  ";NC$(I);n  =  " ; : INPUT ;CI(I,1) 

240  PRINT:NEXT 

250  IF  FL(7)  THEN  260  ELSE  300 

260  INPUT  "Relative  Humidity  (%)  =  ";RH 

270  INPUT  "Input  air  temperature,  C  =  " ;TA:TA=TA+273 

280  P0(7)=(1/760)*EXP(A(7)+BP( 7)/TA) :CO(7)=PO(7)*M(7)*1000000!/(.082055*TA) 
290  Cl (7,1)=C0(7)*(RH/100) 

300  INPUT  "Bed  Temperature, C  =  ",  TE:TE=TE+273:CI(7,1)=CI(7,1)*TA/TE 
310  INPUT  "Bed  Weight,  g  =  ",W 
320  INPUT  "Flow  Rate,  L/min  =  ",F 

330  IF  NFL  THEN  LPRINT  N;"  Component  run  with  ";NS$(FS);"  as  the  sorbent" 
340  IF  NFL  THEN  LPRINT  "Bed  weight  =  ";W;"g"  Flow  rate  = 

";F*1000;"cc/min";"  Temperature  =  ";TE-273;"  C" 

350  IF  NFL  THEN  LPRINT  "Initial  concentrations,  mg/m~3" 

360  FOR  1=1  TO  7:  IF  FL(I)  THEN  370  ELSE  380 
370  IF  NFL  THEN  LPRINT  NC$(I);"  -  " ;CI ( I , 1) 

380  NEXT 

390  FOR  1=1  TO  7:  IF  FL(I)  THEN  400  ELSE  470 

400  P0( I )=(1/760)*EXP(A( I )+BP( I )/TE) :C0( I )=P0( I )*M( I)*1000000! /(.082055*TE) 
410  E( I )=1.987*TE*L0G(C0( I )/CI ( 1,1)) 

420  W0=W0(FS)+L0G(DL( I )/DL( 1) ) 

430  J=l: IF  1=7  THEN  GOSUB  1410 

440  Q( I )=EXP(W0-K(FS)*IE-08*(E( I )/B( I ) )2) 

450  RK(I)=KR(FS)*(H(1)/M(I))-.5:RK(I)=RK(I)*(TE/298)“1.75 
460  Q(I)=1000*(Q(I)/CI(I,1))*(W/F-2.30258*D(FS)/RK(I)) 

470  NEXT 

480  FOR  1=1  TO  7:F0R  J=1  TO  7:W( I , J)=W:NEXT:NEXT 

490  REM  1 Newton-Raphson  Solution  for  liquid  mole  fraction 

500  FOR  J=1  TO  N 

510  FOR  1  =  1  TO  7: IF  FL  ( I )  =  1  then  520  ELSE  530 
520  XC=I 
530  NEXT 

540  X(XC, J)=l: IF  J>1  THEN  550  ELSE  570 

550  FOR  1  =  1  TO  7:  IF  FL  ( I )  THEN  Cl ( I , J)=CI ( I ,  J-l) 
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560  NEXT 

570  FOR  K=1  TO  50 
580  Y=1-X(XC, J) 

590  FOR  L=1  TO  XC-1: IF  FL(L)=1  THEN  600  ELSE  610 

600  Y=Y-(CI(L,J)/CO(L))*(CO(XC)*X(XC, J)/CI (XC,  J) )' (B(L)/B(XC) ) 

610  NEXT  L 
620  DY=-1 

630  FOR  L=1  TO  XC-1: IF  FL(L)=1  THEN  640  ELSE  650 
640 

DY=DY-(CI(L,J)/C0(L))*(C0(XC)/CI(XC,J) )*(B(L)/B(XC))*(C0(XC)*X(XC, J)/CI (XC, J) )- 
“((B(L)/B(XC))-1) 

650  NEXT 

660  DX=Y/DY:X(XC,J)=X(XC,J)-DX: IF  X(XC,J)<0  THEN  X(XC,J)=0 

670  IF  ABS(DX)<. 000001  THEN  740 

680  FOR  L=1  TO  XC 

690  IF  FL(L)=1  THEN  700  ELSE  710 

700  X(L,J)=(CI(L,J)/CO(L))*(CO(XC)*X(XC, J)/CI(XC, J) )~ (B(L)/B(XC) ) 

710  NEXT  L 

720  IF  J>1  THEN  60SUB  1240 
730  NEXT  K 

740  YT=0:IF  NFL  THEN  750  ELSE  760 
750  LPRINT:LPRINT  "Zone  ";J 

760  PRINT.-PRINT  "Zone  ";J  : IF  NFL  THEN  770  ELSE  790 

770  FOR  1=1  TO  7 : IF  FL( I)=l  THEN  LPRINT  NCS(I);"  XL  =  ";X(I,J)  ELSE 

780  NEXT 

790  FOR  1=1  TO  7: IF  FL(I)=1  THEN  PRINT  NC$(I);n  XL  =  ";X(I,J)  ELSE  800 
800  NEXT 

810  FOR  1=1  TO  7:Y(I)=X(I,J)*M(I)/DL(I) 

820  YT=YT+Y(I):NEXT 

830  FOR  1=1  TO  7:Y(I)=Y(I)/YT:NEXT 

840  FOR  1=1  TO  7:  IF  (FL(I)=1)  AND  (X(I,J)>0)  THEN  850  ELSE  930 
850  E(I)=1.987*TE*L0G(C0(I)*X(I ,J)/CI(I , J)) 

860  H0=H0(FS)+L0G(0L(I)/DL(1) ) : IF  Y(I)>0  THEN  W0=W0+L0G(Y(I)) 

870  IF  1=7  THEN  GOSUB  1410 

880  WE(I,J)=EXP(W0-K(FS)*1E-08*(E(I)/B(I))'‘2) 

890  RK(I)=KR(FS)*(M(l)/M(I))'.5:RK(I)=RK(I)*(TE/298)n.75 
900  T(I,J)=1000*(WE(I,J)/CI(I,J))*(W(I, J)/F-2.30258*D(FS)/RK(I) ) 

910  FOR  JJ=1  TO  J-1:T(I,J)=(1/T(I,J)  +1/T(TFL(JJ),JJ))~-1:NEXT 
920  F$*"#####.#  ":G$="#.##" 

930  NEXT 

940  TF=0:F0R  1=1  TO  7 

950  IF  FL(I)=1  THEN  960  ELSE  980 

960  IF  T(I,J)>TF  THEN  TF=T(I,J) 

970  IF  TF«T(I,J)  THEN  TFL(J)=1 
980  NEXT 

990  PFL=0:L=J-1:KF=0 
1000  FOR  K=1  TO  L 

1010  IF  (T(TFL(J),J)>T(TFL(X),K))  THEN  1020  ELSE  1030 
1020  KF«K 
1030  NEXT 

1040  IF  KF  THEN  1050  ELSE  1080 
1050  IF  NFL  THEN  1060  ELSE  1070 
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1060  IF  NFL  THEN  LPRINT  NC$(TFL(J) ) ;"  co-elutes  with  zone  M;KF;"  C  = 
";CI(TFL(J),J);"WE  *  ";WE(TFL(J) ,J) ;"  t/tO  =  ";:IF  NFL  THEN  LPRINT  USING 
G$;T(TFL(J),J)/Q(TFL(J)):PFL=1 

1070  PRINT  NCS(TFL(J));"  co-elutes  with  zone  ";KF;"  C  = 

";CI(TFL(J),J);"  WE  =  " ;WE(TFL(J) , J) t/tO  =  ";:PRINT  USING 
G$;T(TFL(J) ,J)/Q(TFL(J)) :PFL=1 
1080  IF  NFL  THEN  1090  ELSE  1100 

1090  LPRINT,  NC$(TFL( J) ) ;  "  C  =  ";:LPRINT  USING  F$;CI (TFL(J) ,J) LPRINT 
"WE  =  LPRINT  USING  F$;WE(TFL(J) , J) LPRINT  "t  =  LPRINT  USING 
F$;T(TFL( J) ,J) ; : LPRINT  "  t/tO  =  ";:LPRINT  USING  G$;T(TFL( J) , J)/Q(TFL(J) ) 
1100  PRINT, NC$(TFL(J));"  C  =  ";:PRINT  USING  F$;CI(TFL(J) ,J);: PRINT 
"WE  =  ";:PRINT  USING  F$;WE(TFL(J) ,J) PRINT  "t  -  " ; tPRINT  USING 
F$;T(TFL(J) ,J) PRINT  "  t/tO  =  ";:PRINT  USING  GS;T(TFL( J) , J)/Q(TFL(J) ) 

1110  IF  J>1  THEN  GOSUB  1510 
1120  FL(FTL(J))=0 
1130  NEXT  J 

1140  IF  NFL  THEN  LPRINT  CHR$(12) 

1150  END 

1160  DATA  18.6,-3630,18.3,-4180,18.0,-4020,18.6,-4560,18.7,-5030,20.5 
2,-5174 

1170  DATA  6.26,5.03,6.05,5.02,5.98,5.60,6.36,5.92,6.08,5.39 

1 180  DATA  622,2. 14e04 ,421,1. 86e04 ,592,1. 23e04 ,238,1. 7e04 ,534,1. 2e04 

1190  DATA  1,1.25,1.29,1.53,1.76,1.76,0.267 

1200  DATA  96.94,131.39,78.12,92.15,106.17,106.17,18 

1210  DATA  1.218, 1.464, .879, .867, .867, .880, 1.00 

1220  DATA  "1,1-Dlchloroethylene", "Trichloroethylene",  "Benzene",  Toluene", 
"Ethylbenzene",  "o-Xylene" ."Water" 

1230  DATA  "Carbosleve  S-1I" ,"CT" ,"SK-4  V'Carbosleve  G","SpherocarbM 
1240  FOR  1=1  TO  7:Y(I)=X(I,J)*M(I)/0L(I) 

1250  YT=YT+Y(I) :NEXT 

1260  FOR  1=1  TO  7:Y(I)=Y(I) /YT : NEXT 

1270  FOR  1=1  TO  7:  IF  (FL(I)-l)  THEN  1280  ELSE  1390 

1280  E(I)=1.987*TE*L0G(C0(I)*X(I,J)/CI(I,J)) 

1290  W0=W0(FS)+L0G(DL( I)/DL(1) ) : IF  Y(I)>0  THEN  WO=WO+LOG(Y(I)) 

1300  IF  1=7  THEN  GOSUB  1410 

1310  WE(I,J)=EXP(W0-K(FS)*1E-08*(E(I)/B(I))"2) 

1320  RK(I)»KR(FS)*(M(1)/M(I))~.5:RK(I)-RK(I)*(TE/298)~1.75 
1330  T(I ,J)=1000*(WE ( I , J)/CI ( I ,J) )*(W( I , J)/F-2.30258*D(FS)/RK( I) ) 

1340  T(I,J)=(1/T(I,J)  +1/T(TFL(J-1),J-1))"-1 

1350  DC=(2*D(FS)/(2000-D(FS)))*(WE(I,J)-WE(I,J-l))*(T(I,J))/  ( (T (TFL( J- 
1),J-1))-T(I,J)> 

1360  IF  ( DC <0)  AND  (ABS(DC) >CI ( 1,1))  THEN  DC=DC/2  ELSE  1380 

1370  IF  DC<0  THEN  1360 

1380  CI(I,J)  =CI(I,J-1)+DC 

1390  NEXT 

1400  RETURN 

1410  ON  FS  GOTO  1420,1430,1440,1450,1460 
1420  RH=.2:WA=2.46:BA1=.601:BA2=.146:G0T0  1470 
1430  RH= . 365 : WA=2 . 26 : BA1 = . 2163 : BA2= . 0698 : GOTO  1470 
1440  RH=.3:WA=1.59:BA1=.569:BA2=. 1404:G0T0  1470 
1450  RH=.2:WA=2.91:BA1=.724:BA2=. 1437: GOTO  1470 
1460  RH=.38:WA=4! : BA1=. 678: BA2=. 0405 :GOTO  1470 
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1470  IF  CI(7,J)/CO(7)<RH  THEN  1480  ELSE  1490 
1480  W0=W0-WA:B(7)=BA1 
1490  IF  CI(7,J)/CO(7)>RH  THEN  B(7)=BA2 
1500  RETURN 

1510  DC  =  1000*WE(TFL(J) ,J)*(W/F)*(T (TFL(J-l) , J-l)-T (TFL(J) , J) )/ 
(T(TFL(J),J)*T(TFL(1) ,1)) 

1520  C  =  CI(TFL(J),J)  +  DC  :PRINT  "Average  eluent  concentration  = 
";:PRINT  USING  F$;C 

1525  IF  NFL=1  THEN  LPRINT  "Average  eluent  concentration  =  ";:IF  N F L= 1 
1530  RETURN 
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TABLE  B-l.  SUMMARY  OF  SORBENT  PERFORMANCE 

(1,1-Dichloroethylene  at  300  mg/m^  and  5%  RH) 


Sorbent 

Capacity,  mg/ga 

tio.  bb 

sc 

SK-4 

78 

6.85 

14 

Carbosieve  S-ll 

123 

12.87 

5.8 

Graphpac  GB 

0 

(1  min) 

NA 

Carbosphere 

84 

6.47 

17 

Carbopack  C 

0 

(30  sec) 

NA 

WV-B 

32 

1.98 

15 

Spherocarb 

78 

6.45 

12 

LCL 

65 

5.4 

40 

CT 

94 

5.42 

24 

Carbosieve  G 

94 

3.55 

10 

Capacity  based  on  50  percent  breakthrough  time  and  Input 
concentration  level.  (300  mg/ttH  x  flow  rate  x  tsg/bed  weight.) 


bTime  to  10  percent  breakthrough  (note  that  bed  weights  vary). 

cStandard  deviation  expressed  as  a  percentage  of  the  50  percent 
breakthrough  time,  assuming  normal  distribution.  ((t5o-tio)/l*28  x  t5Q.) 


TABLE  B-2.  SUMMARY  OF  SORBENT  PERFORMANCE 

(1,1-Dichloroethylene  at  300  mg/m^  and  85%  RH) 


Sorbent 

Capacity,  mg/ga 

tio.  bb 

sC 

SK-4 

68 

5.63 

13 

Carbosieve  S-ll 

117 

12.3 

6.3 

Graphpac  GB 

0 

(2  min) 

NA 

Carbosphere 

82 

6.77 

12 

Carbopack  C 

0 

(1  min) 

NA 

WV-B 

44 

2.72 

13 

Spherocarb 

99 

7.99 

14 

LCL 

53 

4.38 

4.1 

CT 

110 

7.36 

10 

Carbosieve  G 

88 

3.52 

8.2 

aCapac1ty  based  on  50  percent  breakthrough  time  and  input 
concentration  level.  (300  mg/m^  x  flow  rate  x  tso/bed  weight.) 


bTime  to  10  percent  breakthrough  (note  that  bed  weights  vary). 

cStandard  deviation  expressed  as  a  percentage  of  the  50  percent 
breakthrough  time,  assuming  normal  distribution.  ((t5o-tig)/1.28  x  tsg.) 


TABLE  B-3.  SUMMARY  OF  SORBENT  PERFORMANCE 
(Benzene  at  300  mg/m3  and  5%  RH) 


Sorbent 


Capacity,  mg/g 


Carbosleve  S-ll 
Graphpac  GB 
Carbosphere 
Carbopack  C 
WV-B 

Spherocarb 

LCL 

CT 

Carbosieve  G 


(13  min) 


aCapacity  based  on  50  percent  breakthrough  time  and  input 
concentration  level.  (300  mg/m3  x  flow  rate  x  tso/bed  weight.) 

bTime  to  10  percent  breakthrough  (note  that  bed  weights  vary). 

cStandard  deviation  expressed  as  a  percentage  of  the  50  percent 
breakthrough  time,  assuming  normal  distribution.  ( (t5Q-tio)/1.28  x  t5g). 


TABLE  B-4.  SUMMARY  OF  SORBENT  PERFORMANCE 

(Benzene  at  300  mg/m3  and  85%  RH) 


Sorbent 


K-4 

Carbosieve  S-ll 
Graphpac  GB 
Carbosphere 
Carbopack  C 
WV-B 

Spherocarb 

LCL 

CT 

Carbosieve  G 


Capacity,  mg/ga 

t 10.  bb 

$c 

139 

12.3 

10 

222 

25.7 

4.8 

0 

(10  min) 

NA 

155 

15.1 

2.8 

0 

(3  min) 

NA 

107 

8.1 

1.8 

156 

13.5 

15 

138 

162 

193 

10.6 

11.6 

9.08 

2.5 

3.5 
2.2 

aCapacity  based  on  50  percent  breakthrough  time  and  input 
concentration  level.  (300  mg/m3  x  flow  rate  x  tso/bed  weight.) 

bT1me  to  10  percent  breakthrough  (note  that  bed  weights  vary). 

cStandard  deviation  expressed  as  a  percentage  of  the  50  percent 
breakthrough  time,  assuming  normal  distribution.  ((t5Q-tio)/l-28  x  tso). 
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Figure  B-3.  Water  Adsorption  Isotherm  for  Spherocarb 
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SUMMARY  OF  DATA  FOR  1,1-DICHLOROETHYLENE  (CONTINUED) 
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TABLE  D-l.  SUMMARY  OF  SINGLE-COMPONENT  DATA 
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TABLE  0-1.  SUMMARY  OF  SINGLE-COMPONENT  DATA  (CONCLUDED) 


